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Abstract 

Bose-Einstein condensation has been realized in dilute atomic vapors. This achievement 
has generated immerse interest in this field. Presented is a review of recent theoretical re- 
search into the properties of trapped dilute-gas Bose-Einstein condensates. Among them, 
stability of Bose-Einstein condensates confined in traps is mainly discussed. Static prop- 
erties of the ground state are investigated by use of the variational method. The anlysis 
is extended to the stability of two-component condensates. Time-development of the con- 
densate is well-described by the Gross-Pitaevskii equation which is known in nonlinear 
physics as the nonlinear Schrodinger equation. For the case that the inter-atomic potential 
is effectively attractive, a singularity of the solution emerges in a finite time. This phe- 
nomenon which we call collapse explains the upper bound for the number of atoms in such 
condensates under traps. 



1 



Contents 



1 Introduction 



2 Ginzburg-Pitaevskii-Gross Equation 



2.1 Time-independent case 



2.2 Time-dependent case 



3 Static Properties of Bose-Einstein Condensates 



3.1 Ground state of a Bose-Einstein condensate under axially symmetric mag- 



netic trap 



p. 1.1 Formulation 



p. 1.2 Repulsive case 



p. 1.3 Attractive case 



3.2 



Stability of a two-component Bosc- Einstein condensate 
p. 2.1 Formulation 



P.2.2 Ni = N 2 case 



p. 2. 3 Phase diagrams] 



p. 2. 4 Phase separation 




3.3 Long-range interacting bosons confined in traps 




p. 3.1 Stability of the g 


round state 




p. 3. 2 Charged bosons confined in ion traps 


3.4 Summary 



4 Dynamical Properties of Bose-Einstein Condensates 



4.1 Stability of the D-dimensional nonlinear Schrodinger equation under con- 



fined potential 



4.1.1 Repulsive case 



4.1.2 Attractive case 



4.2 Collapse of the condensate 



4.2.1 4 lme-evolution of the condensate 



4.2.2 Application to 'Li system! 



O Summary 



7 

7 
11 

13 

14 
14 
16 
20 
24 
24 
25 
26 
32 
33 
34 
35 
36 

39 

39 
40 
41 
48 
48 
52 
54 



5 Summary and Concluding Remarks 



56 



2 



A Bose-Einstein Condensation of a Free Boson Gas Confined in Harmonic 



Potentials 



B Pseudopotential 



58 
61 
65 
68 



C The Ground State Energy under the Thomas-Fermi Approximations 
Bibliography 



3 



Chapter 1 
Introduction 



The Bose-Einstein condensation (BEC) 10-01 is one of the most remarkable phenomena 
in quantum many-body systems. The condensation is a logical consequence of the Bose- 
Einstein statistics: a gas of non-interacting bosonic atoms, below a certain temperature, 
suddenly develops a macroscopic population in the lowest energy quantum mechanical 
state. However, BEC occurs at such low temperature that the dimensionless phase-space 
density, p ps = n(Adb) 3 , is larger than 2.612, where n is the number density of atoms, and 
Adb is the thermal de Broglie wavelength, Adb = h/ 1 {2-Kmk-oT) 1 / 2 , with h Planck's constant, 
m the mass of the atom, ks Boltzmann's constant and T the temperature, respectively 
(see, for instance, || as a text book). Accordingly, for many years, BEC had been regarded 
as a mathematical artifact: In December 1924, Einstein wrote to P. Ehrenfest, 'From a 
certain temperature on, the molecules "condense" without attractive forces, that is, they 
accumulate at zero velocity. The theory is pretty, but is there also some truth to it?' §. 
Keesom observed in 1928 the phase transition between He-I (the nomal fluid phase) and He- 
ll (the superfluid phase) in liquid helium, 4 He 0. F. London interpreted this phenomenon 
as BEC in 1938 |J. In liquid helium, however, the nature of the BEC transition is difficult 
to observe clearly because of strong inter-atomic interactions. There exists a non-negligible 
difference between the observed transition temperature T\ = 2.18K and the estimated 
transition temperature T c = 3.14K by an ideal gas approximation. In 1993, the evidence 
of BEC in a gas of excitons in a semi-conductor host has been observed ||. Although 
the interactions in these systems are weak, little information about them is known, and 
thus it is difficult to understand the net effect of BEC in the exciton gas. In 1970's, efforts 
started to observe BEC in a dilute vapor of spin-polarized hydrogen [||, whose inter-atomic 
interaction is sufficiently weak and well understood. However, BEC in the hydrogen gas 
was not achieved for a long time, because of the existence of inelastic inter-atomic collisions, 
which cause trap loss and heating. 

Recently, combinations of various new technologies developed in atomic physics such 
as laser cooling and trapping [10], evaporation cooling and magnetic trap have made it 
possible to increase the phase-space density of a vapor of alkali atoms, and to observe 
BEC of such atoms in controllable situations. In 1995, BEC was observed in a series 



of experiments on vapors of 87 Rb m, 7 Li [IJ, Hj and i3 Na in which the atoms 
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were confined in magnetic traps and cooled down to extremely low temperature of the 
order of 10~ 6 ~ 10~ 7 K (Table [LI] ) [ll|-[]nj. Rigorously speaking, for an ideal gas of 



Table 1.1: Experimental values of parameters, a; s-wave scattering length, T c ; critical 
temperature, N t ; the total number of atoms in trap, u; trap frequency. Note that the total 
number of the atoms is the sum of normal and condensed ones. The values of T c , N t and 
lo are taken from the reference in the right most column. 

iV bosonic atoms in a harmonic potential, the condition n(Adb) 3 > 2.612 is replaced by 
N (huj) 3 j (ksT) 3 > 1.202 where Co = (o^o^u^) 1 ' 3 is the geometric mean of the harmonic 
trap frequencies (Appendix ^). In 1998, BEC of spin-polarized hydrogen atoms was finally 
observed ]20|]. Now, more than twenty groups have succeeded in observing BEC pl]-p4] . 



The experimental realization of BEC of alkali-atoms has stimulated experimental and 
theoretical research in their physical properties. One of the advantages using confined 
atomic vapors is the various choices of atomic systems. The low energy properties are 
characterized by the s-wave scattering length a. While 87 Rb, 23 Na and 1 H atoms have 
positive s-wave scattering lengths [25]-[28], it is known that 7 Li and 85 Rb atoms have 



negative s-wave scattering lengths, corresponding to low-energy attractive interactions [25 
|29l 



In the homogeneous case, it was predicted that the condensate is unstable in three- 
dimension when the inter-particle interactions are attractive P0| . The stability of the 
condensate under magnetic trap has been studied numerically and analytically pT|-^l 



The condensate was predicted to be (meta)stable under the magnetic trap, only when the 
number of atoms is below some critical number. The estimated critical number for 7 Li 
atoms is about 10 3 , which agrees with the observed upper bound of the number of atoms 
in the condensate fT5 |. 

The overlapping condensates of the two different spin states of 87 Rb in a magnetic 



trap were observed |52]|- ||54j| , which has also stimulated theoretical research of the multi- 
component condensates f5^-f6(|. Because of the repulsive interaction between the conden- 
sates and the gravitational force, two condensates were observed to be separated. In the 
study of dynamics of component separation of this system |54]], it was observed that the 
time-evolution is rather complex; the initial configuration quickly damps out and tends to 
preserve the total density profile. With respect to damping, effects of the excitations are 
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not small and it is not yet certain what mechanism is responsible. The time-independence 
of the total density profile is due to the similarity in intraspecies and interspecies scattering 
lentghs in 87 Rb. 

Further, optical trapping of the condensates of 23 Na in several hyperfine states was 
demonstrated |JT|. In a magnetic trap, spins of atoms are polarized. Thus, Bose-Einstein 
condensates are described in terms of scalar wavefunctions (order parameters). In an 
optical trap, on the other hand, degrees of freedom of atomic spin are recovered, and 
consequently the condensate wavefunction behaves as spinor. By applying weak magnetic 
fields with spatial gradient, the degeneracy of the internal degree of freedom is lifted, and 
the freedom of spin orientation gives the formation of spin domains , Q . The domains 
can also be easily miscible by turn off the magnetic fields, which is unlike the experiments 
of 87 Rb j52|. For this system, several theoretical works have also been done 
where Skirmion [63|, spin texture | 65fl , the dynamics of spin mixing 

|68| are studied. 



56 1, spin domain 



formation 



and internal vortex structure 



Keeping the developments in mind, we study static and dynamical properties of the 
Bose-Einstein condensates of atomic gases confined in traps. In Chap. 0, a system of bosons 
with inter-particle interactions of finite range, trapped in an external potential and at very 
low temperature, are considered. Applying the pseudopotential method || |73|] and the 
mean field approach to the system, the Ginzburg-Pitaevskii-Gross equation [p5[|-[f72| is de- 
rived. The time-dependent case, which we call the Gross-Pitaevskii equation, is essentially 
equivalent to the nonlinear Schrodinger equation. In Chap. Bl by use of the variational 



approach |32|, |55fl, static properties of the Bose-Einstein condensates are discussed. The 



ground state properties of a condensate with repulsive or attractive inter-atomic interac- 
tion confined in axially symmetric magnetic trap are investigated (Sec. |3.1| ). The analysis 
is extended to that of the stability of a two-component Bose-Einstein condensate under 
magnetic traps, where the possibility of phase separations is also discussed (Sec. |H7^ ). 
While the inter-particle interaction is assumed to be of finite range in the first two sections 
in this chapter, the properties of the Bose-Einstein condensate of long-ranged interacting 
bosons confined in traps are considered in (and only in) Sec. [3.3| . In Chap. |4], dynamical 
properties of Bose-Einstein condensates are considered. A key idea is the extension of the 
Zakharov's theory [41|. The stability of the wavefunction of the D-dimensional nonlinear 
Schrodinger equation with harmonic potential terms is analyzed (Sec. [4.1| ) and is applied to 
investigate the instability of condensates with effectively attractive inter-atomic interaction 
(Sec. |4.2|) . There, a formula for the critical number of atoms, above which the collapse of 
the condensate occurs, is derived. The last chapter is devoted to summary and concluding 
remarks. 
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Chapter 2 

Ginzburg-Pitaevskii-Gross Equation 



In this chapter, we present the basic equation to describe Bose-Einstein condensates. In 
Sec. |2.1|, following Ref. W% based on the Hartree approximation, we derive the Ginzburg- 



Pitaevskii-Gross equation |p9||-|[72|, which is essentially equivalent to the time- independent 
nonlinear Schrodinger equation. In Sec. |2.2| , we consider the time-dependent case. To 
derive the equation of motion for the condensate, we begin with the second-quantized for- 
mulation. By employing the mean field approach, we obtain the so-called Gross-Pitaevskii 
equation |70|-||72| or equivalently the (time-dependent) nonlinear Schrodinger equation. 



2.1 Time-independent case 



We consider N identical bosonic particles with inter-particle interactions of finite range, 
trapped in an external potential, V(r). We assume that the gas is sufficiently dilute and 
at very low temperature. In such a situation, two-body interaction is dominant and the 
s-wave part plays a central role. We may replace the scattering from an inter-particle 
potential of finite range by that from a hard-sphere of diameter a, which is identical to 
s-wave scattering length in this case. Then, the Hamiltonian operator of the system of 
hard-spheres can be given in certain approximations, the detail of which we explain in the 
Appendix [B], by the pseudopotential Hamiltonian operator J5], 73 1, 

d , 



H=YA-^ i + V(r 




dri 



where m is the mass of a bosonic particle, and 

d 2 d 2 



A, = 

U = 



+ 



+ 



d 2 



dx 2 dy 2 dz\ 



2 ; 



4%h 2 a/m. 



The magnetic trap is well approximated by a harmonic potential, 

m 
~2 



t /" / \ fit f 2 2 i 2 2 i 22 

V(r) = — [u x x + u y y + u z z 



(2.1.1) 

(2.1.2) 

(2.1.3) 
(2.1.4) 

(2.1.5) 
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with (u x , u y , u z ) being trap frequencies. 

In the ground state of the system, almost all bosons may occupy the lowest single- 
particle state because of sufficiently weak inter-particle interactions. Thus, following the 
Hartree approximation, we write the ground state wavefunction, $o( r i; r 2i ••*, r N), i n terms 
of the product of N single-particle state wavef unctions, g(r), 



N 



$o(r*i,T* 2 ,- • -,r N ) = Y[g(vi 



i=l 



where g(r) is normalized as 

(g\g) = Jdr\g(r)\ 2 = l, 
and thus the norm of the wavefunction <3> , defined by (<E>o|$o); 



($o|$o) = J dr i ■■ ■ J dr N \$ (r 



1 • • • r N 



J dr\g(r) 



N 



(2.1.6) 



(2.1.7) 



(2.1.8) 



is equal to unity. From Eqs. Q2.1.1] ) and (|2.1.6| ), we have 
($o\H\Q ) 

= g/ d ^--7 d ^("^ Ai+y(n) ) $c 



lyE/^i-/ dr N %S( ri -r 

i+3 



) — 



(^j$o) 



n 2 



Njdrg*(r) ^-—A + V(r))g(r) 
N(N-l) 



U J dn J dr 2 ^(r 2 )^(ri)<y(ri - r 2 ) [r l2 g(r x )g{r 2 )\ , (2.1.9) 



where the superscript * means the complex conjugate and 

A d 2 d 2 d 2 

A = — + — + 



dx 2 dy 2 dz 2 



;2.1.10) 



In Eq. ( 2.1.9 ), if the product of two single-state wavefunctions g{ri)g{r 2 ) is not singular 
at ?~i2 = 0, the operator (d/dr^) (r i2 ■) can be set equal to unity, which gives 

{%\H\%) = N J drg*(r) (~A + V(r)) g{r) 

N{N ~ %o JdrJ dr'g*(r')g*(r)6(r - r') g{r)g{r>) 



" 2 
N I dr 



9 * {r) ("2^ A + V{r) ) 9{r) + ^ (r)|4 



.(2.1.11) 
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We minimize the functional (&o\H\&o) ( 2.1.11 ) under the constraint: 

($o|$o) = 1- (2-1.12) 

To find the constrained extremum of ($>o\H\$>o), we set the variation of a functional 
($o|-£f|$o) - A^ol^o) equal to zero, 

5 



5g*(r] 



(($o|tf|$o>-/^$o|$o» = 0, 



(2.1.13) 



where \i is a Lagrange's multiplier. Substituting Eqs. Q2.1.8|) and (|2.1.11| ) into Eq. fl2.1.13|) 
we obtain 

- ^A^(r) + 7(r)^(r) + (iV - l)C/ |<?(r-)| 2 ^) = m {v). 
2m 



Introducing a wavefunction ^(r) 



(r) = iV 1/2 #(r 



we get 



with 



2m 



V{r)^(r) + 1 



N 



U |#(r)| 2 #(r) = (j&(r) 



(2.1.14) 
(2.1.15) 



= / dr|^(r)| 2 = N. 



(2.1.16) 
(2.1.17) 



The wavefunction \1/ has many names; an order parameter, a macroscopic wavefunction 
and a Bose-Einstein condensate wavefunction. For sufficiently large N, we can neglect the 
1 / iV-order term appearing in the left hand side of Eq. fl2.1.16|) . The result is 

h 2 

A* (r) + V(rU(r) + U \^(r)\ 2 *(r) = utf (r). (2.1.18) 

2m 



We refer to Eq. ( |2.1.18| ) as the Ginzburg-Pitaevskii-Gross equation [B9"|-||72|| with an ex- 
ternal potential term. This equation may also be referred to as the time- independent 
nonlinear Schrodinger equation. 

From Eqs. Q2.1.11 ) and ( |2.1.15| ), the ground state energy of the system, E, is 

r *2 



E 



N / dr 



|- \Vg(r)\ 2 + V(r) \g(r)\ 2 + (iV - 1) ^ \g(r) 





~h 2 


1 dr 


2m 



|W(r)| 2 + V(r) |*(r) 



1 - 



2 

l_\ U 

N 



;2.i.i9) 



which we call the Ginzburg-Pitaevskii-Gross energy functional f6"9]|— 172| with an external 
potential term. On the other hand, multiplying the both sides of Eq. (|2.1.16j ) by ty*(r) 
and integrating it, we have 

r *2 



Jdr |-|W| 2 + ^)|*| 2 + (l- 1)^ |¥| 



E+{1--)^ /dr|* 



(2.1.20) 
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Another derivation of Eq. ( [2.1.20]) may be instructive. We give a small change to the 
single-particle wavefunction g(r), and denote the resultant by g(r), 

g(r) = (l + e)^g(r), (2.1.21) 

where e is a sufficiently small constant and g(r) satisfies Eqs. ( |2.1.7|) and ( 2. 1. 14j) . Similar 
to ( |2.1.21|) , we also define 

V(r) = (1 + e) 1/2 *(r), (2.1.22) 
where * satisfies Eq. ( |2.1.16| ). From the normalization condition of \&(r) ([2.1.17 ), we have 

(*|*)= [ dr\^(r)\ 2 = N + eN. 



(2.1.23) 

In ( 2.1.23|) , we can regard e N as a virtual change of the number of particles, 5N, 

5N = eN. (2.1.24) 
From (2.1.6) with fl2.1.21| ), we construct a transformed wavefunction of the system, $ , 



N 



$o = \{~g{ri), 



(2.1.25) 



Then, the norm of the wavefunction $o is 

($o|$o> = 



dr\g(r) 



N 



:i+e 



ss 1 + eN =1 + 5N. (2.1.26) 

From ( |2.1.12| ) and ( ^.1.26| ), we obtain a variation of the norm of $o> 

8(%\%) = (%\%) - ($ |$o) = SN. (2.1.27) 

The corresponding change of the energy is calculated as follows. By replacing \I/ with \I> in 
( g 1.19|) , we get 



/ 



dr 



2m 



V{r) 



! at; 2 



E + 5E, 



where 



| dr [^|V»|' + V(r)|*|"+(l-l)l7.|«| 



'2.1.28) 



(2.1.29) 



Because * satisfies Eq. ( |2.1.16|) , the variations 5N (|2.1.27|) and 5£ (|2.1.29| ) satisfy the 
extremum condition ( |2.1.13|) , and thus we have 

5E-fi5N = 0. (2.1.30) 

Substituting ( f2.1.29|) into ( |2.1.3q ), we get 



5E _ 1 

In ~n 



dr 



EL 

2m 



|W| 2 + ^(r) |*| 



N 



u m 



(2.1.31) 



which is the same as ( [2.1. 20| ). From ( [2.1.31[ ), it is clear that /i, introduced as a Lagrange's 
multiplier, has a meaning of the chemical potential of the system. 
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2.2 Time-dependent case 



In this sub-section, we consider the time-dependent case, \I> = \fr(r,t). To derive the 
equation of motion of ty, we use the second-quantized formulation. Let ^f(r, t) and W (r, t) 
denote the bosonic annihilation and creation operators, respectively. The (equal-time) 
commutation relations among the operators are 



4f(r,t)4?(r r ,t) - $f(r',t)i?(r,t) = 0, 
6(r - r'), 



and the second-quantized Hamiltonian of the system, H, can be written as 



H = dr 



h 



Ua 



V*t . + V(r)&$ + ^t^t^ 
2m 2 



(2.2.1) 
(2.2.2) 



(2.2.3) 



The time-evolution of the operator &(r, t) obeys the Heisenberg equation, 



•« 9 f 
*** 



(2.2.4) 



Substituting the Hamiltonian (|2.2.3| ) into Eq. Q2.2.4 ) and using the commutation rela- 
tions (ggj ) and (gjj) , we get 



<5 . ft 2 . » ... 

ift— tf = + V(r)ty + U &^. (2.2.5) 

<9t 2m v 7 v 7 

We denote an expectation value by ( • ). The Heisenberg equation for the bosonic 
operator fl2.2.5|) gives 







(*) = -^A (*) + V(r) (*) + U ($tM) . 



(2.2.6) 



According to the mean field theory [15, TBI , we may replace expectation values of the 
bosonic annihilation and creation operators by the condensate wavefunction ty(r,t) and 
its complex conjugate ^*(r,t) respectively, 

(*(r, t)) = *(r, f), (*t(r, *)) = **(r, *). (2.2.7) 

For the third term in the right hand side of Eq. ( |2.2.6 ), we take the following approximation, 

(^M) « (\j/ f ) (£\ = |^(r,t)| 2 ^(r,t). (2.2.8) 

Substituting ( gg77p and (|2"X8"1) into Eq. flETBD , we have 



i^-^(r, t) = -^A*(r, *) + V(r)^(r, t) + U |tt(r, t) | 2 (r, *), (2.2.9) 
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which is called the Gross-Pitaevskii equation ]7D1- ||72|| with an external potential or the 
(time-dependent) nonlinear Schrodinger equation. The one-dimensional nonlinear Schrodinger 
equation is known to be integrable and has been studied extensively related to various areas 
of physics 75 1. The equation ( |2.2.9|) can be written in a variational form, 



in—— 

dt 



5 



5^* 



Em, 



;2.2.10) 



where the functional E[-] is defined by Eq. (|2.1.19| ) with the 1/iV-order term deleted. We 
note that, by setting 

tf(r,t) = exp(-i/xt) W(r), (2.2.11) 

in Eq. (|2.2.9|) , we obtain the Ginzburg-Pitaevskii-Gross equation (|2.1.18 ) again. 

The existence of non-zero *&(r, t) = (^(r, t)^ has an important meaning in physics, the 
breakdown of the gauge symmetry. One of the consequences is that the condensate has 
a definite phase 9(r,t) as defined by ty(r,t) = f(r, t) exp[i 9(r, t)] where / and 9 are real 
functions. 
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Chapter 3 



Static Properties of Bose-Einstein 
Condensates 



In this chapter, we consider the static properties of Bose-Einstein condensates. In Sec. 3.1 



by using the variational approach, we examine the ground state properties of a Bose- 
Einstein condensate with repulsive or attractive inter-atomic interaction confined in axially 
symmetric magnetic trap Employing a gaussian trial wavefunction to describe the 

condensate, we derive the minimum conditions of the Ginzburg-Pitaevskii-Gross energy 
functional ( |2.1.19p . In the repulsive inter-atomic interaction case, it is shown that, if the 
trap has a high asymmetry and thus the equi-potential surface of the trap is "cigar" or 
"pancake" shaped, we can assume different approximations to each of axial and transverse 
directions. We also compare the energy with the one obtained from the Thomas-Fermi 



approximation p3| . In the attractive interaction case, different from the repulsive case, 
the energy function has only a local minimum. The local minimum of the energy no longer 
exists above some critical number of atoms. This can be regarded as a static theory of the 
collapse. We calculate the critical number for three shapes of the equi-potential surface of 
the trap: almost sphere, cigar and pancake. 

In Sec. |3.2| , we extend the analysis in Sec. [JTT]to discuss the stability of a two-component 
Bose-Einstein condensate under magnetic traps |K| . We first investigate the case where the 



numbers of two species of atoms are the same. The effect of interspecies interaction on the 
critical number of atoms is explicitly shown. We also present phase diagrams for various 
combinations of the scattering lengths. We further develop the variational approach and 
consider the phase separation of the condensate. 

Inter-particle interactions of finite range are well approximated by the delta function. 
We may think of the long-range interaction like a Coulomb potential, because of the possi- 
bility of Bose-Einstein condensation for all bosonic atoms and molecules with and without 
charges in future. Then, in Sec. |3.3| (and only in this section), we consider the properties 
of the Bose-Einstein condensate of long-ranged interacting bosons confined in traps. As 
in Sec. |3.1| , we derive the minimum conditions of the energy, and obtain the ground state 
energy. It is interesting to see that the condensate of long-range interacting bosons under 
traps is stable both for repulsive and attractive cases. 
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3.1 Ground state of a Bose-Einstein condensate under 
axially symmetric magnetic trap 

The ground state properties of a Bose-Einstein condensate with repulsive or attractive 
inter-atomic interaction confined in axially symmetric magnetic trap are studied through 
the variational approach f4"4" |. It is shown that, if the trap has a high asymmetry and thus 



the equi-potential surface of the trap is of "cigar" or "pancake" shape, the properties of 
the condensate are drastically changed from the case of isotropic (spherically symmetric) 
trap. 

3.1.1 Formulation 

We consider the ground state properties of a Bose-Einstein condensate confined in axially 
symmetric magnetic trap, 

y(r) = ^iri + c^ 2 ). (3.1.1) 

Here, the axis of the symmetry is chosen to be the z-axis, and r_|_ denotes the radius of the 
projection of the position vector r on the xy-plane, m the atomic mass, and u z and u± 
the trap (angular) frequencies along the z-axis and in the xy-plane, respectively. 

Following Baym and Pethick |33f and Fetter |J2| for isotropic case, we employ as the 



ground state wavefunction a gaussian trial function, 



N 



1/2 



*(r±,*) = ^-^[4 J «p[-(ri/4+^K)/2], (3.1.2) 

where d± and d z are variational parameters. The macroscopic wavefunction of the conden- 
sate, \l/(r), is assumed to be axially symmetric, as the magnetic trap is. This assumption 
is reasonable as far as the system is near the ground state. We have normalized the trial 
function (|3.1.2|) to the number of the particles in the condensate, N, 



/*oo /*oo 

2vr / r ± dr ± / dz\^{r ± , z)\ 2 = N. (3.1.3) 

JO J-oo 

For the weakly interacting system at zero temperature, the number of the non-condensate 
particles is considered to be negligible. The parameters d± and d z measure the extent 
of the wavefunction in the radial and axial directions. If there is no interaction between 
atoms, the exact ground state wavefunction is obtained by setting d± and d z equal to 
characteristic oscillator lenghs, l± and l z respectively, which are defined by 

i± = (— ) 1/2 . ^m 1/2 . (3^» 



In the Hartree approximation (see Sec. |2.1|) , the ground state energy of the system is given 



by a Ginzburg-Pitaevskii-Gross energy functional |69|-|72l with the harmonic potential 
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terms, 







1 dr 


— |W| 




2m 



m 



uj*rl + uiz z )m' + 



2nh 2 a, 



1^1 



(3.1.5) 



where a is the s-wave scattering length. We note that the minimization of the functional 
E — fiN, where /i is the chemical potential, gives the Ginzburg-Pitaevskii-Gross equa- 
tion |6{|-|72] with the harmonic potential terms, 



Substituing ( |3.1.2|) into (|3.1.5|) , we obtain 

-2 



Anh 2 a 



m 



\^| 2 \P = ;U\I>. 



E(d_i_,d z 



N 



Tr 



2mdj 



m 
~2 



+ '^u 2 d\ 



+ 



N f h m , 
7 ( + 2^*, 



h 2 aN 2 



(2Ti) l / 2 m d\d z 



(3.1.6) 



(3.1.7) 



The kinetic energy and the potential energy are respectively proportional to d 2 and d 2 . 
We introduce dimensionless variational parameters 



s± = l±/d±, s z = l z /d z , 
and dimensionless experimental constants 

5 = ouJujx, G ± = (2/7r) 1 / 2 Na/l ± . 



(3.1.8) 



(3.1.9) 

The constant 5 measures the anisotropy of the trap. The constant G± represents the 
strength of the interaction and is positive (negative) for the repulsive (attractive) inter- 
action case. We may use other constants such as G z = {2/n:) 1 / 2 Na/l z ^ but those can be 
expressed in terms of 5 and G±. Substituting ( |3.1.8|) and ( |3.1.9|) into (|3.1.7| ) yields 



E(s±, s z 



Nhu± 



[sl + s-^ + ^sl + s^ + S^G^ 



(3.1.10) 



We note that the large (small) s corresponds to the contraction (spread) of the condensate. 
We investigate the existence of a minimum of E(sj_, s z ) as a function of the two variational 
parameters s± and s z . The minimum conditions, 



dE 

ds± 
dE 

ds 7 . 



0. 
0. 



(3.1.11) 
(3.1.12) 



determine the location of the minimum. Substituting (gTlCD into (|3.1.11|) , we get 



K s ± - s ± 



^ + 5 1/2 G ±S± s z = 0. 
Similarly, substitution of Q3.1.10D into ( |3.1.12|) gives 



(3.1.13) 



s z - s. 



5- l l 2 G { s 2 =0. 



(3.1.14) 
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We also calculate the second derivatives and d f^g s of E(s±, s z ), 

= Nhuj ± (l + 3sl A + 5 l/2 G ± s z ), 

US j_ 

= ^Nhuj ± (l + 3s:% 



o 2 e 



5 1/2 Nhuj ± G ± s ±J 



ds±ds z 

from which we have the determinant of the Hessian matrix, A, 
d 2 E d 2 E 



A 



d 2 E \ 2 



ds\ ds 2 z \ds±ds z J 

5(Nhuj ± ) 2 Ul + 3sl 4 + 5 1/2 G ± s z )(l + 3s; 4 )/2 - G\s 2 _ 



(3.1.15) 
(3.1.16) 
(3.1.17) 



(3.1.18) 



In the following, we deal with the minimum conditions ( |3.1.13| ) and fl3.1.14|) analytically 
through appropriate approximations depending on the values of S and G±. 



3.1.2 Repulsive case 

We consider the case where the inter-atomic interaction is repulsive. In this case, the s- 
wave scattering length a and accordingly the constant G± defined by ( |3.1.9| ) are positive. 
From (|3.1.13| ), we have 



- 1 



-S 1/2 G i si s 7 < 0. 



(3.1.19) 



Then, the parameter s± is smaller than 1 at the minimum point. In a similar manner, we 
can show from ( |3.1.14| ) that s z is smaller than 1 at the minimum point. Those indicate 
that a repulsive inter-atomic interaction always acts to expand the condensate size fl32| . 
From ( |3.1.15| ), we see that |^ is positive, and by substituting ( |3 . 1 . 1 3| ) and ( |3.1.14| ) into 
( P.1.181 ), we have 

A = 5(Nhcu ± ) 2 (3sl 4 + s~ A + 5s7>; 4 - 1), (3.1.20) 



which gives 



A > 5{Nhu ± ) 2 {sl 



I) > 0. 



(3.1.21) 



Here, we have used the fact that s± and s z are smaller than 1. Therefore, we confirm that 
a point (sj_, s z ) which satisfies ( |3.1.13|) and (|3.1.14j ) is indeed a minimum point. 
1 ) Weak interaction approximation 

We first consider the case that the third terms in both (|3.1.13|) and (|3 . 1 . 1 4 ) are negli- 
gible, which means the following conditions, 



S_L - S 



> 8 1/2 G ±S±S ; 



sz 3 » d-^G^s 2 



(3.1.22) 
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This approximation corresponds to the situation where the inter-atomic interaction is suf- 
ficiently weak in both xy- and z-directions. In this case, the approximate solutions of 
(gTTD and (j3XTp are 

a± = 1 - ^ 1/2 G ± , a x = l- ^<T 1/2 G ± . (3.1.23) 

Recall that s± = s z = 1 recovers the exact result for non-interaction case. By substituting 
( 3.1.23 ) into ( 3.1.22 ), we obtain a criterion of the approximation for 5 and G±, 

G± <min{<T 1/2 ,5 1/2 }. (3.1.24) 

We refer to (|3.1.24j) as the weak interaction case (region I in Fig. |3.1.1| ). From (|3.1.8 ), we 



Figure 3.1.1: Approximations for the repulsive inter-atomic interaction: I. weak interaction 
case (|3.1.24) i II. strong interaction case (|3.1.31| ), III. intermediate case-1 (|3.1.40|) , IV. 
intermediate case-2 ( 3.1.52 ). The abscissa and ordinate represent the anisotropy of the 
trap, 5, and the strength of the interaction, G±, respectively. 



have the original variational parameters as 



1 + h^Gj] l ± , d z =(l + -S-^G,] I 



d i 



which give the aspect ratio of the condensate, A = d z /d±, 



A = 5~ x ' 2 



l + ^r 1 / 2 -^/ 2 )^^ 



(3.1.25) 



(3.1.26) 



We note that, if 8 < 1 {5 > 1), the aspect ratio ( |3.1.26|) is larger (smaller) than that for 
the non-interaction case, G± = 0. The energy E of the condensate (|3.1.10|) is estimated as 



E 



Nhuj± 



2 + 5 + 5 1/2 G A 



(3.1.27) 



2) Strong interaction approximation 
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We consider the case that the first terms in both (p.1.13 ) and ( |3 . 1 . 14| ) are small. The 
conditions are 



s± < S 1/2 G±s ± s z ~ s 



3 , s z « S-^G , s 2 ~ s -3. 



_L ) 



(3.1.28) 



This is similar to the Thomas-Fermi approximation |33| where the effect of the inter-atomic 
interaction is assumed to be dominant over the kinetic one. In this case, the approximate 
solutions of s± and s z are given by 



5 1/2 (5G ± Y 



1 + -(S 2 - 3)(5G A 



1/5 



1 



5 V 2 



-4/5 



(SG A 



-4/5 



(3.1.29) 
(3.1.30) 



By substituting (|3.1.29| ) and (|3.1.30| ) into (|3.1.28| ), we obtain a criterion of the approxima- 
tion for S and G±, 

G ± ^>max{5-\5 3/2 }. (3.1.31) 

We call fl3.1.31|) the strong interaction case (region II in Fig. |3.1.1| ). The corresponding 
values of d±, d z , A and E are calculated as 



d ± 
d z 
A 
E 



(SG, 



,1/5 



-4/5 



/ 



-4/5 



v-1 



6 



l + -(5 2 -l)(5G ± )- 4 / 5 



(3.1.32) 
(3.1.33) 
(3.1.34) 



(SG ± ) 



2/5 



1 + - 1 + 



(SG A 



57T 



3/5 



j2/5 ^vriVaX 2 / 5 



-4/5 



(*G_l) 



-4/5 



(3.1.35) 



The energy (|3.1.35| ) is equal to Eq. (7) in Ref. |33| to the leading order. It is interesting to 
compare the leading order term of the energy (|3.1.35|) with the energy obtained from the 
three-dimensional Thomas- Fermi approximation (see Appendix 0), 



E, 



3D 



5 /15g\ 
7 V8vr; 



2/5 



muz 



2\ 1/5 



muj. 



2\ 1/5 

y \ 



into 



2\ 1/5 



(3.1.36) 



We see that both of which are proportional to iV 7 / 5 . 
3) Mixed approximations 

In both weak and strong interaction cases, we have treated similarly the two minimum 
conditions, Q3.1.13| ) and ( |3.1.14| ). However, there exist situations where we can give different 
approximations to each condition. 

First, we consider the case that the third term in ( |3.1.13| ) and the first term in fl3.1.14| ) 
are small, satisfying 



s± ~ si 3 > 5 1/2 G ± s ± s z , s z < 5- 1/2 G ± s 2 ± ~ sf. 



(3.1.37) 
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This implies that the effect of the inter-atomic interaction is negligible in ^-direction, but 
important in z-direction like the Thomas-Fermi approximation. In this case, s± and s z are 
approximated as 



{5- 1 ' 2 G_ 



1-1/3 



3 o 



(3.1.38) 
(3.1.39) 



By use of ([3.1 .38Q and (|3.1.39| ) in (|3.1.37 ), a criterion of the approximation is given by 

S 1/2 < G ± < <T\ (3.1.40) 

which requires 

<J«1. (3.1.41) 

Therefore, this approximation is allowed only when the equi-potential surface is elongated 
in the ^-direction and thus "cigar" shaped. We call (|3.1.40| ) the intermediate case-1 (region 
III in Fig. |3.1.1[ ). Correspondingly, the original parameters d± and d z , the aspect ratio A 
and the energy E are calculated as 



d± 
d z 
A 
E 



1 + - 4 (5G ± )^ 



(d^G^ 3 [l + \5- l l 2 G L )-^ - -(5G ± ) 2 / 3 
3 6 



r l/2 (r l/2 G±) l/3 N + l (r l/2 G±) -4/3 _ l ((5G±) 2/3 

3 X 2 



> i, 



Nhu ± 



2 + ^G ± ) 2 / 3 



The energy can be rewritten as 



3 • 2 1 / 3 

E = Nhuj± + — {mu ± g/h) 2/3 

OTT 



2\ 1/3 



JV5/3 



where 



(3.1.42) 
(3.1.43) 
(3.1.44) 
(3.1.45) 

(3.1.46) 
(3.1.47) 



g = 4irh 2 a/m. 

The first term is the ground state energy of two-dimensional harmonic oscillators. We note 
that the second term is proportional to iV 5 / 3 , in the same way as the energy obtained from 
the one-dimensional Thomas- Fermi approximation (see Appendix |C|), 



E 



ID 



2/3 



2\ i/3 



N 5/3 . 



(3.1.48) 



Second, we consider the opposite case; the first term in (|3.1.13|) and the third term in 
( p.!.14[ ) can be treated perturbatively, which means the following conditions, 



s ± < 5 1/2 G±s±s z - s 



3 

_L ) 



sz 3 » d-^G^s 2 



(3.1.49) 
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In this case, the effect of the inter-atomic interaction is dominant in ^-direction but 
negligible in z-direction. The parameters s±_ and s z are given by 



16 4 



(3.1.50) 
(3.1.51) 



Substitution of ( |3.1.50|) and Q3.1.51]) into ( |3.1.49|) gives a criterion of the approximation for 
S and G± as 

<T 1/2 < G± < <5 3/2 , (3.1.52) 

which requires 

<J»1. (3.1.53) 

Therefore, this approximation is allowed only when the equi-potential surface of the trap 
is an extremely flat spheroid, like a "pancake" . The condition (|3.1.52j ) is referred to as the 
intermediate case-2 (region IV in Fig. |3.1.1| ). The expressions of d±, d z , A and E are given 
by 



d± 

d z 
A 
E 



(5 l/2 G±) l/4 L _ _L (r 3/2 G±) l/2 + h S l/* G A-l 

16 4 



l + I( r 3/2 G±) l/2 



S- 1 l\8 l ' 2 G A .)- 1 ' A \l + —(S-^GA 1 / 2 - -(S^GA- 1 

16 4 

Nhuj? 



<1, 



2 

NTiuj, 



\ + 2{5- z ' 2 G 1 _) 1 " 



+ 



muj 



2 \ 1/2 



2 1 /4 7I -3/4 



1/4 



1/2^3/2 



(3.1.54) 
(3.1.55) 
(3.1.56) 

(3.1.57) 



with g defined by ( |3.1.47| ). We note that the second term in the last expression of ( |3.1.57| ) 
is proportional to N 3 ^ 2 , just like the two-dimensional Thomas- Fermi approximation energy 
(see Appendix 0), 



2 /2\ 1 /2 



2D 



3 Vtt 



2\ 1/4 



2\ 1/4 

y \ 



1/2^3/2 



(3.1.58) 



while the first term is the ground state energy of one-dimensional harmonic oscillators. 



3.1.3 Attractive case 



We investigate the case that the inter-atomic interaction is effectively attractive; the s-wave 
scattering length a and accordingly the constant G± defined by ( |3.1.9| ) are assumed to be 
negative. There, the energy function (|3.1.10| ) has only a local minimum, different from the 

Since G 



previous section 32 



is negative, we have from ( 3.1.13 ) 

A -1) = _5l/2 G|S 4 > 0, 



(3.1.59) 
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which shows that s± is larger than 1 at the minimum point. Similarly, s z is larger than 1 
at the minimum point. We thus observe that an attractive inter-atomic interaction works 
to shrink the condensate size ||32| . 

By use of the minimum condition ( 3.1.13] ) in ( 3.1.15Q , we find that |^ is always positive 
at a local minimum point, 



d 2 E 
dsj 



ANTiuj\s~i > 0. 



(3.1.60) 



If G± is sufficiently small, the minimization of the energy ( |3.1.10|) is attained at (s±, s z ) ~ 
(1, 1), and we see that, by substituting s± ~ 1 and s z ~ 1 into (|3.1.20|) , the determinant of 
the Hessian matrix A is positive at this point. However, it can be shown that A becomes 
zero for some value of G±, at which the local minimum of the energy disappears, leading 
to the instability of the condensate. This instability will be discussed in Chap. | as a 
dynamical problem. The critical value of G± is determined by taking the value of A at the 
minimum point to be zero, 



A 



<5(iV7^) 2 (3sJ_ 4 + s~ A + 5sJ>; 4 - 1) 



0. 



(3.1.61) 



with the minimum conditions ( gTTl and (|TT|). From (glial) , (MM) and ( HH) , 
we get a set of three equations which determine the stability condition, 



:4+5)(4-3)(4-i) 

a si +5 



2„4 



(3.1.62) 
(3.1.63) 



G ± = -5- 1/2 s-\l - si 4 ). (3.1.64) 

We can solve ( ]3.1.62|) analytically for s 4 . However, the expression is too complicated 
to get useful information. In the following, we investigate Eqs. Q3.1.62| )- ( p. 1.64] ), by using 
some approximations to them. 
1 ) Almost spherical trap 

Here we consider the case that the trap is almost spherically symmetric, namely 



5 = lo z /uj± = 1 + e, |e| <g 1. 
In this case, Eq. ( |3.1.62j ) and then Eq. (|3.1.63|) are approximately solved as 

2e\ 



- - 51/1 + 1) 



5 1 / 4 1 - 



- . 



9 J 



From (3.1.66), we get the expressions of the original variational parameters, 



5- 1 / 4 1 



4 = 5-^(i + |) l Z) 



which give the aspect ratio A = d z jd\ 



A 



1-1. 
6 



(3.1.65) 



(3.1.66) 



(3.1.67) 



(3.1.68) 
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Using ( |3.1.66| ) in ( |3.1.64| ), we obtain the critical value of G±, denoted by G_l )C , 



G 



5 5 /4 



e 



(3.1.69) 



The result agrees with the one obtained by Fetter (Eq. (7) in Ref. |[32||), to the leading order. 
Further, by substituting (|3.1.9|) into (|3.1.69|) , we have the critical number of particles, N c , 
at which the condensate becomes unstable, 



Nr. 



2 3 / 2 tt 1 / 2 h 
5 5 / 4 la 



1 - 



6 



0.671- 



23/2^1/2 ^2/3^/3 

5 5 /4 | a | 



(3.1.70) 



We also calculate the (locally) minimized energy at G± = Gj_ c , defined as E c , 



5 l/2 



Nhui 



5 l/2 



(3.1.71) 



In an experiment for 7 Li [|T^], the values of frequencies u z and u± were ui z /2tt pa 117 Hz 
and u±/2tc pa 163 Hz. The s-wave scattering length of 7 Li was observed to be a = — 27.3ao 
where a-o is the Bohr radius |p9| |. Using the data in the formula (|3.1.70| ), we have A^ c = 1450. 
This value agrees with the experimental result, iV c = 650 ~ 1300 [13]. 

2) Cigar shaped trap 

We consider the case that the equi-potential surface of the magnetic trap is "cigar" 
shaped, 

5 = e< 1. 



In this case, Eqs. ( |3.1.62|) and ( |3.1.63| ) are approximately solved to give 



s^^{l+ e -\, s z = 
From (|3.1.73| ), we get d±, d z and A as follows, 
d± = 3~ 1/4 (l - jj l ± , d z 



4\ 1/4 



-i/ 2 ! + 



3e2 
16 



(3.1.72) 



(3.1.73) 



.4 



3\ x / 2 
2, 



1 + 



gy/4 

16 j ' 



(3.1.74) 



(3.1.75) 



We observe that <i z is about the order of l±_ rather than l z , and accordingly the aspect 
ratio of the condensate A is about the unity. It is intriguing that, near G± = G± tC , the 
condensate shape does not reflect the anisotropy of the trap, which is assumed to be highly 
elongated in ^-direction. 

From ( |3.1.64| ) and ( |3.1.73| ), we obtain G_l jC , 



G 



±,c 



2 l/2 



i + 



ljk 2 

16 



(3.1.76) 
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and then the critical number of particles iV c , 



71 



1/2 



3 3/4 



1 + 



13e 2 \ / 



13e 2 \ Z, 



w- it" - 778 1 + w" e- 



(3.1.77) 



This critical number is larger than the one obtained in the almost spherical case ( |3.1.70|) , 
if we use the same l± and a in both cases. The critical (locally) minimized energy E c is 
obtained as 



Er 



3I/2 



1 - ^) Ntos 



(3.1.78) 



3) Pancake shaped trap 

We consider the opposite case to the previous one; the equi-potential surface of the 
trap is an extremely flat spheroid, like a "pancake" , namely 

e 1 » 1 



5 = e - x >l (e<l). 
In this case, the approximate solutions of Eqs. ( |3.1.62| ) and ( |3.1.63| ) are 

e 2/3 



From (13.1.800 , we get 



8 ± = 4 1 '*e- 1 '*, s z = l + 



4 -l/4 e l/6 Z = 4 -l/4 e -l/3, , 



(3.1.79) 



(3.1.80) 



(3.1.81) 



(3.1.82) 
V z / 

Note that d± is the order of e 1 / 6 /^ or equivalently e~ 1//3 / z , meaning that d±_ is much larger 
than l z but much smaller than Zj_. On the other hand, d z remains the order of l z . Therefore, 
the aspect ratio A = d z /dj_, which is the order of e 1//3 , is much larger than the ratio between 
the characteristic oscillator lengths, l z /lj_ = e 1 / 2 . We observe again that the condensate 
shape near the instability is not highly anisotropic. 
Substitution of (|3.1.8C| ) into (|3.1.64j ) gives 



from which we get 



Nr. 



71 



1/2 



(V,, = - e V2fl_^2/3 



2 l/2 



1_ V/ 3 )^r« 1.25(1 -V/M^r 



4 



4 



(3.1.83) 



(3.1.84) 



We see that the critical number ( |3.1.84|) is larger than the ones obtained in the previous 
cases (|3.1.7q) and ( |3.1.77| ), if we set l z in ( p.l.84Q equal to l± in (|3.1.7q) and ( |3.1.77| ), for 
a fixed a. The critical minimized energy E c is calculated as 



Er 



1 + e 4/3 )iVftu; 2 



(3.1.85) 
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3.2 Stability of a two-component Bose-Einstein con- 
densate 



The stability of a two-component Bose-Einstein condensate is an interesting subject ex- 
perimentally and theoretically Using the variational method, we investigate the stable 
and unstable conditions as functions of the particle numbers and the s-wave scattering 
lengths |60|. We find that interaction between the different species of atoms (simply, in- 
terspecies interaction) has significant effects on the stability of each condensate. We also 
consider the phase separation of the condensate. 



3.2.1 Formulation 

We consider a two-component Bose-Einstein condensate in isotropic harmonic trap poten- 
tials. We denote by m i; cjj, N iy and a„ the mass, the trap frequency, the number of atoms, 
and the s-wave scattering length of the z-th component respectively. The energy functional 
of the macroscopic wavefunctions ^(t*), which are normalized to iVj, can be written as 



£[*i,* 2 ] 



j dr 






i=l,2 \ 



2m,- 



\V%(r)\ 2 + Vi(r)\^i(r) 



+ 



v-* 27ih 2 au 

lv,:(r] 



i=l,2 



rrii 



4+ 2^ 
m 12 



\Mr) 



where VAr) 



1, 2) stand for the isotropic potentials, 

1 



V,r) 



2 2 

-vriiUJ-r . 
2 1 



(3.2.1) 



(3.2.2) 



Minimizing the energy functional flOp with (ggg) gives the Ginzburg-Pitaevskii-Gross 
equation for binary mixtures of the condensates under magnetic traps. The last term in 
( p.2.1|) represents the interactions between the species- 1 and the species-2 atoms, whose 
strength is proportional to the s-wave scattering length a 12 between them and is anti- 
proportional to the reduced mass m,\ 2 = m\m 2 j{m\ + m 2 ). The positive (negative) value 
of dij means the effectively repulsive (attractive) interaction. To perform variational cal- 
culations, we use trial functions of the gaussian form 

*i(r) = N?*-*dk* exp (~ j . (3.2.3) 



The parameter di measures the extent of spreading of ^(r). By substituting ( fj.2.3|) into 
( p.2.1[ ), one gets the energy as a function of variational parameters d\ and d 2 

h 2 

12) X X ' 




+ ^27r^a« (27r) _3 iV 3 d _3 



+ 



2-Kh a 12 3 AT AT / j9 , 9 
m 12 v 



(3.2.4) 
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To investigate the effects of all the inter-atomic interactions as rigorous as possible, we con- 
centrate on a case where mi = m 2 = m and uji = uj 2 = lu. We introduce the dimensionless 
energy e, scattering length a^, and width parameter Aj as 

E = lL' aij = ~P Xi = 7' (3.2.5) 



where I = Jh/(mw) (the characteristic oscillator length). Using e, a^, and Aj instead of 
E, aij, and dj, we rewrite ( pj.2.4 ) as 



e = ?EW + A*) + 4=E^««V 3 + -^iV!iV 2 a 12 (A? + A, 2 )" 1 - (3.2.6) 



With respect to the global properties of the function s = e(Xi, A 2 ), we remark the following: 

(i) A global maximum does not exist because e — > oo as Ai, A 2 — > oo. 

(ii) The existence of a minimum depends on the signs of the scattering lengths a^. 
The energy e = e(Xi, A 2 ) takes a local minimum when Ai and A 2 satisfy 

0C : 1 - Af 4 - \ -N iail \f ~ ^N 2 a 12 (Xj + A 2 )" 1 = 0, (3.2.7) 



d\i V Ti" \A" 



cfe / 2 8 

: 1 - A 2 ~ 4 - J-N 2 a 22 X 2 5 - _JViai 2 (A? + A 2 )^ = 0, (3.2.8) 



d\ 2 z U 0F 

and the eigenvalues of the Hessian matrix are positive: 

d 2 e d 2 e ( d 2 e \ 2 „ d 2 e 



Therefore we can examine the existence of a stable state with gaussian wavefunctions using 
( |3.2.7| )-( |3.2.9|) . We analyze these conditions for various physical situations described by 



and iVj, as follows. 



3.2.2 Ni = N 2 case 



We begin with the case where N\ = N 2 = N and an = a 22 = a. We have set an = a 22 to 
express the results in analytic forms. The conditions (|3.2.7|) and ( |3.2.8|) give 



Ai = A 2 , Xl- Ax = J-N(a + a l2 ). (3.2.10) 

V vr 



Figure |3.2.1| illustrates the stable region for a, ai 2 and N, where (Ai, A 2 ) satisfy ( |3.2.9|) and 



( p.2.10| ). The unstable region implies the instability (collapse) of the condensates. From 
the analysis of the boundary curve, we obtain a critical number N c (= N± c = N 2c ), which 
gives an upper limit to the number of atoms, and the parameter A c (= Ai c = A 2c ), 



N c = - 2 3 1 , A c = 5~* (a 12 <0, a + a 12 <0), (3.2.11) 

52 a + a i2 
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Figure 3.2.1: Stability of the two-component condensate: Ni = N 2 case for (a) a > 0, (b) 
a < 0. The boundaries of stable and unstable regions (the solid lines) are determined by 
Eqs. Q3.2.9D and ( |3.2.10| ). The dashed lines are discussed in |3.2.4 
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-i 



1 + 4^ 



A, 



1 + 4^) (ai2>0, a<0). (3.2.12) 
\a\ 



The conditions (|3.2.11|) and ( 3.2.12 ) are worth to be discussed. First, if there is no inter- 
action between two species of atoms, that is, a 12 = 0, we recover N c = 2\[2tx j^\a\ for 
the one-component case |32|]. Recall that the isotropic limit of ( |3.1.70|) gives this N c . Sec- 
ond, when the interspecies interaction is effectively attractive (repulsive), that is, otyi < 
(«i2 > 0), the critical number of atoms N c is reduced (enhanced). Thirdly, it is remarkable 
that the instability occurs even for a > if ai 2 < —a. 



3.2.3 Phase diagrams 

We construct phase diagrams for N\ and iV 2 : find the boundary line N 2 = f(Nx) which di- 
vides iViiVVplane into stable and unstable regions. In other words, we pursue the maximum 
of N 2 function of N±. 

1) an > 0, «22 > 

This case corresponds to a situation that both intraspecies interactions are repulsive. If 
interspecies interaction is not strongly attractive, the condensate is stable for any Ni and 
-/V 2 because of their intraspecies repulsions. Therefore we consider here the case ai 2 < 
and | cti2 1 is sufficiently large. We also consider the case that {Ni,N 2 }, {Ai,A 2 }, and 
{«n,a 22 , |«i 2 |} are sets of the same order quantities. We represent as N, A and a the 
quantities which are the same order as {Ni, N 2 }, {Ai, A 2 } and {an, a 22 , |ai 2 |} respectively. 
We start from the following assumption for (Ni,N 2 ) on the boundary N 2 = f(Nx) and 
other variables, and we check the validity of the assumption at the end. 
The assumption: on the boundary, f(N\)/Ni remains finite as Ni tends to infinity, and Ai 
and A 2 do not get too large when Ni, N 2 3> 1, so as to satisfy 

1 w 0{\- 4 ) < 0{NaX~ 5 ). (3.2.13) 
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This assumption is reasonable since we expect that the attractive interspecies interaction 
suppresses the growth of A as iV gets large. 

Under the assumption ( [3.2. 13|) , neglecting the first and second terms in fl3.2.7|) and 
( |3.2.8| ), we obtain the ratios N 2 /Ni and A 2 /Ai in the case a\ 2 < — y/oina 2 2, 



N2 
iVi 



7 



|Ol2| 
«22 



1 ± 



\ 



«11«22 
, l«12| 2 . 



(3.2.14) 



F12I 

aiiQ!22, 



1 ± 



V 



\ 



Q11Q22 

, l«12| 2 . 



(3.2.15) 



/ 



Equation (3.2.14) determines asymptotic lines of the boundary curve for sufficiently large 
Ni and N 2 . Here we remark that de/dXi = do not give a useful information about the 
boundary line (or the critical number), but with the approximation used here, de/dXi = 
imply A = 0. That is the reason why the ratio f(N\)/Ni has been pursued in ( |3.2.14 ). 
To calculate Ai (not A2/A1), we must take account of terms in the first order which are 
neglected here. By substituting A2/A1 = {3(1 + 6) and N 2 /Ni = 7(1 + 5), where \S\, |e| <C 1, 
into (|3X7|) , ( gXg) and A = 0, we obtain 



Ax = 5 



_i / a n {3 3 + a 22 ^ 
\au(3 7 + a 2 27, 



(3.2.16) 



This shows that Ai is a quantity of 0(1) and therefore the assumption ( [3.2.13D is confirmed. 

Next we consider the case of small N\. For simplicity, we consider the case where 
«u = a 2 2 = ex. The boundary curve N 2 = f(Ni) passes through (Ni,N 2 ) = (N C ,N C ), 
where iV c = 2y / 27r/54 (|a 12 | — a) ( |3.2.11| ). We examine the behavior of the boundary curve 
in a neighborhood of this point. We substitute Ni = N c + SNi and Aj = A c + 5Xi into 
de/dXi = and A = 0. The zeroth order terms give A c = 5~4 and the first order terms 
give 





«12 








• a 



(SXx - 5X 2 ) 



1 



[a6Ni + a 12 5N 2 ) 



(ai 2 SNt + a5N 2 ) 



Therefore, we get 



dN 2 
diV"! 



-1 at (N 1 ,N 2 ) = (N C ,N C ). 



(3.2.17) 



(3.2.18) 



Figure |3.2.2| illustrates the phase diagram based on the above discussion. We see that the 
large attractive interspecies interaction overcomes the repulsive intraspecies interaction and 
leads to the collapse. We expect the similar tendencies for an 7^ «22- 
2) a u < 0, a 22 < 
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Figure 3.2.2: Phase diagram for an = a 22 = a > and ai 2 < —a. 



This case corresponds to a situation that both intraspecies interactions are attractive. 



The curve N 2 = f{N\) intersects with the iV 2 -axis at (Ni, N 2 ) = (0,N^), where N : 



r(0) 

2c 



2v27r/54 |ct 22 1 is obtained from the one-component theory. We first consider the behavior 
of N 2 = f(Ni) near this point. Setting Aj = limjv^o \ and substituting Aj = A^ + 5Xi 
into ( 3.2.7| ) , ( p . 2 . 8|) and A = 0, one gets from the zeroth order terms 



1-A? } -2 



9 «12 
2 

I "22 I 



Af 



i 

5~*. 



(3.2.19) 
(3.2.20) 



and from the first order terms 
dN 2 



diVi 



n 5 "12 
2a 

"22 



5lAf )2 + l 



1-Ai' 



(o)- 



at (iVi,iV 2 ) = (0,iV 2 ( c 0) ), (3.2.21) 



where Xf 1 ^ is determined by ( |3.2.19|) . 

Next, in the case an = a 22 = a < 0, we consider the neighborhood of the point 
(JVi, N 2 ) = (JV C , N c ), where N c is given by (|3.2.11|) and ( |3.2.12|) . If a 12 < 0, N c is the same 
as the case where a > and ai 2 < —a. Therefore, 



diV 2 
diVi 



-1 at (N 1 ,N 2 ) = (N C ,N C ). 



(3.2.22) 



In the case ai 2 > 0, substituting N{ = N c + 5Ni and Aj = A c + 5Xi into ( |3.2.7|) , (|3.2.8|) and 
A = gives 



"12 

a 



— 5X\ — 5X 2 



1 



2,71 



(a6Ni + a l2 5N 2 ) 



{a 12 5Ni + a5N 2 



(3.2.23) 
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Figure 3.2.3: Phase diagram for a u = a 22 = a < 0. (a) a\ 2 > 0, (b) a < a\ 2 < 0, (c) 
ol\i < a < 0. 



Thus, we get 

diVo 

—^ = 1 at (N l7 N 2 ) = (N c ,N c ). (3.2.24) 
dNi 

Figure |3.2.3| illustrates the phase diagrams for various a\ 2 . We can see that the repulsive 
(attractive) interspecies interaction enlarges (narrows) the stable region. 
3) a u > 0, a 22 < 

In this case, the species-1 atoms have repulsive interactions and the species-2 atoms have 

l 5 

attractive interactions. In a neighborhood of (Ni,N 2 ) = (0, 2v27r/54 |a 22 |), the behavior 
of the boundary N 2 = f(Ni) is the same as that in the previous case. We investigate here 
limjv^oo f(Ni). First we assume N 2 is sufficiently smaller than Ni. Setting N 2 = in 
( gXTp and (|3^8|) gives 

Af - Ai - ^Niau = 0, (3.2.25) 

1 - A 2 - 4 - -^N iai2 {\l + Al)"i = 0. (3.2.26) 
V 71 " 

As seen from the condition ( |3.2.25 ), the width Ai gets larger as Ni increases because of the 
repulsive intraspecies interaction. If a\ 2 < (attractive interspecies interaction), we have 
A 2 < 1 from ( |3.2.26| ) and A 2 Ai is satisfied. This condition is also satisfied for positive 
and small a 12 . When 1 < A 1; A 2 < A 1( and N 2 < N u one gets from ( |3.2.7|) and ( |3.2.8|) 





= 0, (3.2.27) 

V 71 

2 AT 2 a 22 A 2 - 5 - 4=^iai2Ar 5 = 0. (3.2.28) 

7T 



Since 



|| = §JVi ( 1 + 4^iV 1 « 11 Ar 5 j = yiV x > 0, (3.2.29) 
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d 2 F 60 

= Nl N 2 a 12 \~ e \ 2 « 0, 
0X10X2 v 71 " 



and 

we can use d 2 ejdX\ = instead of A = 0. The condition d 2 e/dX 2 = gives 

1 + 3A2 4 + 4 J-N 2 a 22 X 2 5 - -^AWAr 5 = 0. 

V 7T V7T 

From ( |3.2.27| ), ( |3.2.28j ) and ( |3.2.31| ), we obtain the critical number of N 2 , 

1 

' 4 



(3.2.30) 



(3.2.3i; 



N, 



2c 



2\/2k I / 22^12 



5* «22 



an 



(3.2.32) 



This implies that N 2c is independent of jVj. and thus limjVi-^oo N 2c /Ni = 0. 

If «i2 is positive and large enough, we expect that the species-2 condensate broadens 
in its size and N 2c /N\ does not converge to zero as N\ tends to infinity. We can pursue 
N 2c for a fixed large Ni without using A = directly, as follows. We assume X 2 ^> 1 for 
sufficiently large ol\ 2 and obtain from Eqs. ( p.2.7| ) and ( E>.2.8| ) 



2 8 
-N x a lx X{ h + — A^ 2 ai2 (X\ + A 



7T 



_ 5 
2"\ 2 
2 



g 

2a 22 A2" 5 + ^A^icua (x\ + A 
v 71 



_ r, 
2\ 2 
2, 



(3.2.33) 
(3.2.34) 



Equations p.2.33j ) and ( ^.2.34| ) give a relation for N 2 /N\ as a function of the ratio (3 
A2/A1, 

-«11+«12 (^)~ § 



N2 



(3.2.35) 



|«22|/3- 5 + «l 2 (^)^ 

Figure |3.2.4| illustrates a function in ( |3.2.35| ). In the region f3 ~ 0, 7(/3) ~ {(—an + 
22a 12 )/|a 2 2|}/? 5 - Therefore when a 12 > 2^a n , the maximum of 7(/3) determines the 
critical number of N 2 . Explicit form of max^7(/3) can be given in two special cases where 
( = 2^a\ 2 ja\\ is 1+ and 00. 
(i) Setting ( = 1 + e and expanding in e give 



2\ 2 «n 7 
e 2 . 



— = max 7 /5 r_ 
N t P \7J \a 22 \ 

(ii) Supposing ( is large and expanding in l/( give 



(3.2.36) 



N 2c ,m , 1 I , '>22 

— = max 7 /3 ~ 1 - - U + 
iVi /3 C I V «n 



(3.2.37) 
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Figure 3.2.4: Function defined in ( |3.2.35| ). We set a n : a 22 ■ ol 12 = 1 : -4 : 2~2. 



Figure 3.2.5: liniAr^oo f(Ni)/Ni as a function of 22a 12 /an- We set an : a 22 = 1 : —4. 



Figure 3.2.6: Phase diagram for an > and a 22 < 0. (a) a\ 2 > 2 a an, (b) < «i2 < 

_ 5 

2 2 an, (c) «i2 < 0. 
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Summarizing the above discussion, limjVx-KX) f(Ni)/Ni as a function of 22cei 2 /an is illus- 
trated in Fig. |3.2.5 . 

Figure p . 2 . 6| shows phase diagrams for various ai 2 . If a.\ 2 > 2~2atn, N 2c gets large 
proportional to N± by enlarging the width X 2 . On the other hand, in the case a% 2 < 2~2ctn, 
the extent of the species-2 condensate remains small because of the attractive intraspecies 
interaction and its critical number converges to a finite number as N\ increases. 



3.2.4 Phase separation 



In the experiments of two-component BEC |52|- |M , phase separation due to the gravita- 
tion and the repulsive interspecies interaction was observed. In this sub-section, we discuss 
the possibility of phase separation caused by the latter. For simplicity, we consider the 
case where 

Nx = N 2 = N, a u = a 22 = a, oj x = u y = uj±, uj z = 5uj^_. (3.2.38) 

First, we assume that the anisotropic parameter 5 satisfies 8 > 1 (pancake-like trap). 
To discuss a possibility of phase separation, we generalize the trial functions into 



tyi(r) = N"n ±(d x d y d z ) 2 exp 



x 



2dl 2dl 2dl 



P ± V2Pi 



X 



d x J ' 



(3.2.39) 



Here + (— ) corresponds to i = 1 (i = 2), and real parameters Pq and P\ satisfy the relation 
Pq + P\ = 1, which is required by the normalization condition, 



J dr|^i(r)| 2 = N. 



(3.2.40) 



The non-zero value of P\ detects the phase separation. Substitution of the trial functions 
( |3.2.39|) into (|3.2.1| ) gives the energy as a function of the parameters d x , d y , d z and P\. In 
terms of dimensionless variables, 



E _ a a 12 . _ d x _ d y _ d z 

£ — 7 5 a — i—j a i2 — ~r~, *x — -j—, *y — — , \ z — — , 

nuj± l± lj_ l± l± l z 

where l± = \Jhj (muj±), l z = sjft,/ (muj z ), the energy becomes 

N 



(3.2.41) 



£ = 2 



1 + 2Pt){\i + k>) + a; + x- 2 + 5(K + k 

+^K% 1 K 1 S 1/2 {aN (-\pt + Pi + l) 



+ai2 N _ 3 p2 + 1 

The necessary condition for phase separation is 



0, — = at < Pi < 1, A fc > (k = x, y, z). 



dP{ 



(3.2.42) 



(3.2.43) 
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Thus, the conditions: 



de 



0. 



de 



pf=o 



(k = x,y,z), 



(3.2.44) 



p2=o 



jive the boundary curve of phase separation. Solving (|3.2.44j) gives 

/ 2a 12 \* 



A. T — A,, 



\azi2 — a 



^ = S»J-N(avi - a), 



N 



0jr J a + a 12 + \J (a + a 12 ) 2 + 85 2 a 12 (a 12 - a) 



(3.2.45) 



(3.2.46) 



2 S(ot\2 — a)i(2ai 2 ) 4 

If a±2 is larger than a and N is larger than the value (|3.2.46|) , phase separation occurs. 



Second, we consider the case, 5 < 1 (cigar-like trap). Using trial functions 



m i (r) = N^K *(d x d y d z ) * exp 



x 



V 



2dl 2d 2 2^)( Po± ^ Pl i 



(3.2.47) 



we can calculate the critical number of phase separation, 



N 



2tx ^2 (a + a 12 ) + JS 3 (a + a l2 ) 2 + 85ai 2 (ai 2 - a) 



(«12 - a)t(2ai 2 )4 



(3.2.48) 



In the limit 8 — *■ 1, (|3.2.46| ) and ( |3.2.48|) coincide. The boundary curves of the phase 
separation for this case are shown in Fig. |3.2.1| (the dashed lines). We see that the larger 
interspecies interaction causes phase separation for the smaller number of particles. 



3.3 Long-range interacting bosons confined in traps 

Throughout this paper except this section, we consider a system of neutral bosonic atoms 
under magnetic traps. The interaction between neutral atoms is well approximated by 
the delta function, V(r) = g5(r). We may think of the other limit, that is, the long- 
range interaction like a Coulomb potential, V(r) = g/\r\. In principle, we believe that the 
Bose-Einstein condensation occurs for all bosonic atoms and molecules with and without 
charges. 

In this section, we consider the properties of the Bose-Einstein condensate of long- 
ranged interacting bosons confined in traps. The problem is rather academic, but turns 
out to be quite interesting: The condensates in the short-ranged and long-ranged cases 
have different properties. 
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3.3.1 Stability of the ground state 

We again start from the Ginzburg-Pitaevskii-Gross energy functional, 

2m 



dr 



^— |W| 2 + -mwV |*| 2 



+ \J dT I dr ' U ( r - r ') l^( r )! 2 l*( r ') 



where the trap is assumed to be an isotropic harmonic potential, and 

U(r) = gj \r\ . 



(3.3.1) 
(3.3.2) 



The coupling constant g can be positive or negative in this sub-section. 

We assume that the ground state wavefunction \l/(r) depends only on r = \r\. With 
this assumption and an integral formula for a fixed r, 



-i/2 _ J 2/r (for r' < r) 
~~ \ 2/r' (for r' > r), 



(3.3.3) 



sin 6 d6 [r 2 + r' 2 — 2rr' cos 6 ) 

the last term in ( fj.3. 1| ) is written as 

[dr I dr' - — - — - |^(r)| 2 |^(r')| 2 = 16tt 2 ^ f° dr r 2 \^{r)\ 2 f°° dr' r' \^{r')\ 2 . (3.3.4) 
J J \r — r'\ Jo Jr 
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We calculate the ground state energy by the variational method. As a trial function, 
we choose 



^r(r) = Cexp 



2d? r 



(3.3.5) 



where C and d are real constants to be determined. The particle number N and the ground 
state energy are calculated as 



N 
E 



C 2 ^ 2 d\ 

3n 3 / 2 h 2 C 2 , 37T 3 / 2 muj 2 C 2 



Am 



4 2 y 



(3.3.6) 
(3.3.7) 



By use of ( p. 3. 6 ), we eliminate the normalization constant C in ( 3.3.7|) . The result is 



E(d) 



3h 2 N 3mu 2 Nd 2 N 2 g 



+ 



+ 



Amd? ' 4 ' (2tt)W 
It is convenient to introduce a dimensionless parameter A, 

A = d/d , 

where d = [/i/(mw)] 1 '' 2 . Then, ( |3.3.8|) is rewritten as 



E(X) = -Nhu 



(3.3.8) 



(3.3.9) 



(3.3.10) 
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a 



(3.3.1i; 



where 

^2 Ng fmui\ 1 / 2 
n %u \ h J 

We can show that E(X) has an absolute minimum irrespective of the sign of a, that is, 
the condensate of long-range interacting bosons under harmonic traps is stable both for 
repulsive and attractive interactions. 

It is interesting to compare the above result with that for the delta function case. The 
same calculation with U{r) = g 5(r) gives 



E d (\) 



-Nhu 
2 



2 V 



+ A^ 



2 mgN f muj\ 1 / 2 
7T A-k% 2 \~h) 



(3.3.12) 
(3.3.13) 



A clear difference between ( |3.3.10| ) and ( |3.3.12| ) is the A-dependence of the interaction 
term. The powers of A, A -1 in (|3.3.10| ) and A -3 in ( |3.3.12|) , are easily understood from the 
scaling property of the interaction potential U(r) under the transformation r — > Xr. In 
the case of the long-range interaction, the kinetic term proportional to A~ 2 is dominant 
over the interaction term proportional to A -1 as A — > 0. This explains why the collapse of 
the condensate does not occur for the long-range interacting bosons confined in traps. 



3.3.2 Charged bosons confined in ion traps 

We consider the Bose-Einstein condensate of charged bosons confined in ion traps. Based 
on the results obtained in the sub-section [3.3. 1| , we further examine the ground state energy. 
Setting g = e 2 , e being the electric charge, we have 



E(X) 



-Nhcu -(A" 2 + A 2 ) + a c A- 1 



where 



(To 



^ Ne 2 /m^y/ 2 

7T TlUJ \ tl J 



(3.3.14) 



(3.3.15) 



We note that a e is positive. 

We minimize -E"(A) with respect to A. The condition, dE/dX = 0, gives 

3A 4 — 3 — a c X = 0. (3.3.16) 

For a weak interaction case where a e is small, the approximate solution of (|3.3.16|) is 
A « 1 + (0-/12). Using this in ( |3.3.14fl , we obtain 

E 



3 -m, (1 + £ 1 



(3.3.17) 



In the non-interaction limit a e — > 0, the exact result E = 3Nhu)/2 for harmonic oscillators 
is recovered. 



35 



For a strong interaction case, a e ^> 1, ( |3.3.16| ) is approximately solved to have 



A = 3- 1/3 a c 1/3 + a' 1 . (3.3.18) 
Then, the ground state energy is obtained as 

E = ^3 4/3 N?aua 2/3 + ^3 5/3 NhLua^ 2/3 . (3.3.19) 

The dominant term in Q3.3.19| ) shows the ^-dependence of the energy explicitly as 

o4/3o2/3 

E = ±_l I7 _ m l/3 w 2/3 e 4/3 JV 5/3 > ( g 3 2Q) 

This result corresponds to the Thomas-Fermi approximation in the sense that the kinetic 
energy is negligible. There should be no confusion with the formula ( |(J.6| ) where the 
interaction is the delta function type. 

3.4 Summary 

We have investigated the static properties of Bose-Einstein condensates. In Sec. |3TT| , we 
have studied the ground state properties of a Bose-Einstein condensate with repulsive or at- 
tractive inter-atomic interaction confined in axially symmetric magnetic trap. The gaussian 
trial wavefunction has two variational parameters which are determined by minimum con- 
ditions of the Ginzburg-Pitaevskii-Gross energy functional with harmonic potential terms. 



For the repulsive case (sub-section |3.1.2| ), we find four situations according to which terms 
appearing in ( |3.1.13| ) and (|3.1.14 ) can be neglected in a first approximation. The ap- 
proximation conditions in each situation are summarized in Fig. |3.1.1| with the anisotropy 
parameter of the trap 5 = uo z /uj±_ and the relative strength of the inter-atomic interaction 
G± = (2/ir) 1 / 2 Na(muj±/h) 1 / 2 . For all cases, we can estimate the variational parameters 
s_i_ and s z or equivalently dj_ and d z , the aspect ratio of the condensate A = d z /d± and the 
energy E at the minimum point. In the cases 1) and 2), we have used common approx- 
imations to both of the two minimum conditions Q3.1.13| ) and ( |3.1.14| ). For the case that 



the inter-atomic interaction is sufficiently weak (resp. strong) in both xy- and z-directions, 
the approximation condition is given by ( |3.1.24| ) (resp. ( [3.1.31 )) corresponding to the re- 



gion I, the weak interaction case (resp. the region II, the strong interaction case). The 
energy ( |3.1.35| ) and the one obtained from the three-dimensional Thomas-Fermi approxi- 



mation have the same iV-dependence, proportinal to N 7 ^ . It is interesting that, when the 
anisotropy parameter 6 is much larger or smaller than unity, there exist cases where we 
can give different description to each direction. In the case 3), we have first considered 
the case that the effect of the inter-atomic interaction is negligible in ^-direction but 
important in ^-direction. This occurs when the trap is cigar-shaped. The approximation 
condition is obtained as (|3.1.40|) , corresponding to the region III, the intermediate case-1. 
The energy ( |3.1.46| ) consists of two parts. The first term describes the two-dimensional 
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harmonic oscillator. The iV-dependence of the second term is iV 5 / 3 which is the same for 
the one-dimensional Thomas- Fermi approximation. Next, we have considered the case that 
the effect of the inter-atomic interaction is dominant in ^-direction but negligible in z- 
direction. This is realized when the trap is pancake-shaped. We obtain the approximation 
condition as ( 3.1.52Q , corresponding to the region IV, the intermediate case-2. We find 
here that the second term in the energy ( |3.1.57|) is proportional to iV 3//2 , as the energy 
obtained from the two-dimensional Thomas-Fermi approximation. To summarize, in the 
above results, the N- dependence of the energy is consistent with the geometrical shape 
of the condensate and the dimensionality of the Thomas-Fermi approximation. For the 
attractive case (sub-section |3.1.3| ), the energy function ( p.l.lOp has only a local (not an 
absolute) minimum. As we expect, the local minimum of the energy disappears for some 
critical value of G±, which implies the instability (collapse) of the condensate. By setting 
the value of the determinant of the Hessian matrix A at the minimum point to be zero 
( [3. 1.61 ) with the minimum conditions ( p. 1.13 ) and (p. 1.14 ), we calculate the critical value 
G± <c and obtain the critical number of particles N c . First, for the case that the trap is 
nearly spherical, N c is estimated as ( 3.1.70|) , which is essentially the same as the one ob- 
tained by Fetter [32[]. Next, for the cigar-shaped trap, N c , estimated as ( |3.1.77| ), is larger 
than the critical number for the nearly spherical trap. The aspect ratio of the condensate 
A at G± = G± tC ( |3.1.75| ) is the order of unity. Similarly, for the pancake-shaped trap, we 
obtain the critical number as ( |3.1.84|) , which is larger than the ones for the nearly spherical 
and the cigar-shaped traps. 

In Sec. |3.2| , we have investigated the stability of a two-component Bose-Einstein con- 
densate. For various choices of "dimensionless" s-wave scattering lengths Oy, we have con- 
structed the phase diagram in iViiVVplane where Ni is the number of species-z (i = 1,2) 
atoms. Recall that is related to the s-wave scattering length by ( |3.2.5|) . We sum- 
marize the salient points of the results, paying attention to the collapse. 

1) an > 0, a22 > 0. A sufficiently strong attractive interaction between the different species 
can overcome repulsive interactions within each condensate, which leads to their collapse. 

2) an < 0, a22 < 0. A repulsive (attractive) interaction between the different species makes 
the critical number of atoms larger (smaller) than the one for one-component case. 

3) an > 0, a 22 < 0. For a sufficiently strong repulsive interspecies interaction, the criti- 
cal number of species-2 atoms, increases as the number of species-1 atoms, N\. For 
an attractive interspecies interaction, N^c decreases and approaches a finite value as Ni 
increases. 

A basic assumption of our analysis has been the gaussian trial functions ( |3.2.3| ). While 
the gaussian ansatz is reasonable for the coexisting condensates in magnetic traps, there 
is no guarantee that ( |3.2.3| ) exhausts all the possible stable states. Thus, in the sub- 
section |3.2.4| , we have modified the trial functions, and discussed the possibility of the 
phase separations. In the pancake trap case, we generalize the trial functions into ( |3.2.39p , 
and obtain the critical number of atoms ( 3.2.46Q , above which the phase separation occurs. 
In the cigar trap case, the trial functions are modified as (|3.2.47p , and the critical number 
of atoms for the phase separation is ( |3.2.48| ). The boundary curves of the phase separation 
in the isotropic limit are plotted by the dashed lines in Fig. |3.2.1 . 
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In Sec. |3.3| , we have considered the properties of the Bose-Einstein condensate of long- 
ranged interacting bosons confined in traps. As in the other sections in this chapter, 
we have chosen a gaussian function (|3.3.5|) as a trial function, and obtained the energy 
function ( 3.3.1 0| ) . The energy ( 3.3. 10| ) has an absolute minimum irrespective of the sign of 
the strength of the interaction, a, and therefore the condensate of long-range interacting 
bosons under traps is stable. By use of the minimum conditions of the energy ( |3.3.16| ), we 
have investigated the ground state energy. For the weak interation case, the energy ( [3.3.171) 
is approximately equal to that for harmonic oscillators. For the strong interaction case, the 
leading term of the energy is proportional to N 5 ^ 3 . The results obtained in this sub-section 
may give useful information, when the Bose-Einstein condensation for bosonic ions under 
ion traps will be observed in future. 
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Chapter 4 



Dynamical Properties of 
Bose-Einstein Condensates 



In this chapter, we investigate the dynamical properties of Bose-Einstein condensates. In 
Sec. |4.1| , we first investigate the stability of the wavefunction of the D-dimensional non- 
linear Schrodinger equation with harmonic potential terms ||76|| . For the repulsive case, 
we prove that the wavefunction is absolutely stable. For the attractive case, by extending 
the Zakharov's theory JF7, 78|, we show that the singularity of the wavefunction surely 



emerges in a finite time when the total energy of the system is negative, and even when 
the energy is positive, the wavefunction collapses in a finite time for a certain class of the 
initial conditions. In Sec. |4~2| , based on the results obtained in Sec. [O], we investigate the 
dynamics of the condensates analytically, and prove that the singularity of solution emerges 



in a finite time when the total energy of the system is both negative and positive ||41||-||43 
In the analysis, the initial wavefunction is assumed to be gaussian with two parameters. 
We present a formula for the critical number of atoms. Further, by improving our analysis 
quantitatively, we apply it to the assembly of 7 Li atoms. Then, within a reasonable pa- 
rameter region, the estimated critical number for 7 Li atoms is comparable with the upper 
bound for the number of atoms in the recent experiment of BEC. 

4.1 Stability of the D- dimensional nonlinear Schrodinger 
equation under confined potential 

We investigate the stability of the wavefunction of the D-dimensional nonlinear Schrodinger 
equation with harmonic potential terms. The analysis in this section is a mathematical 
introduction to the problem of the collapse of the condensate in Sec. [O. 



Let x = (xi, x 2 , ■ ■ ', Xd) be a point on a D-dimensional Euclidian space R/°. As usual, 
time is denoted by t. Then the equation for the wavefunction t) is expressed as 

2i^ = -A0 + V(x)<f> + g\(j>\% (4.1.1) 
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Here A = Y, n =i d /dx^ is the D-dimensional Laplace operator, g = ±1, and 



D 



V(x) 



E2 2 



V n > 0. 



(4.1.2) 



71=1 



The coupling constant g = +1 (g = —1) represents the repulsive (attractive) self-interaction. 
The numerical factors in ( |4.1.1| ) with ( f4.1.2| ) are chosen for convenience. 

Equation (|4.1.1|) without external potential (V(x) = 0) has been studied by many 
researchers [[77[j -|^6||. Zakharov J77] and Zakharov and Synakh |7H[ studied the multi- 
dimensional nonlinear Schrodinger equation (D = 2,3) for the attractive case (g = — 1) 
as models of Langmuir waves in plasma and light waves in nonlinear media. They showed 
that, when the total energy of the system is negative, the singularity surely emerges in 



the wavefunction in a finite time. The phenomenon is termed the collapse. Weinstein |80 



derived the condition that the collapse occurs even when the energy of the system is 
positive. Then, it is extremely interesting to examine how the condition obtained by 
Weinstein is changed when the harmonic potential terms ( |4. 1 .2| ) are added. 



4.1.1 Repulsive case 

We deal with the repulsive case, namely g = +1, and show that the wavefunction <p(x,t) 
in Eq. ( |4.1.1| ) is absolutely stable for D > 1, when the wavefunction is confined in the 
harmonic potential (|4.1.2j ). 

Equation ( f4.1.1| ) has two constants of motion, 



iV 
H 



2 dx, 



\V<P\ 2 + V{x 



2 + fl0| 4 



dx, 



(4.1.3) 
(4.1.4) 



corresponding to the conservations of particle number N and energy H, respectively. Here 
and hereafter, dx = dxi ■ ■ ■ dxD- We define a norm, II • |L, as 



\F(x)\ p dx 



1/p 



(4.1.5) 



For g = +1, (|4T4D with PX^ ) reads 

H = /(Wi 2 + e 

\ n=l 



vy n \<p\ 2 +^\ 



dx 



D 



n<r\\2 



+ 



(4.1.6) 



n=l 



Because the third term in the second line of Eq. (^4.1.6|) is positive, we have inequalities, 

D 



h > + e ui\\x n <f>wi > wvm + a 2 im|2, 



(4.1.7) 



71=1 
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where the constant A is the smallest among {u n }, 

A = min(K}). (4.1.8) 
Further, by applying the Cauchy-Schwarz inequality (the uncertainty relation in physics), 

IIV0HI \\z4Wl > (^||0||^/2) 2 = (DN/2) 2 , (4.1.9) 

to ( [4.1. 7|) , we obtain 

H > \\V(f)\\l + {kDN/2f\\V(j)\\f > ADN. (4.1.10) 

From ([4.1. 10|), we see that ||V0||2 is bounded, and thus the wavefunction <fi(x,t) is abso- 



lutely stable for the repulsive case. 
4.1.2 Attractive case 

Next, we consider the attractive case, g — — 1. To investigate the time-evolution of the 
wavefunction 0(x, t), we introduce the expectation value of the square of the radius, \x\ 2 = 
x\ + x\ + • • • + x 2 ^, with respect to </>, 

(\x\ 2 (t)} = \\x(f)(t)\\l = J \x\ 2 \(f)(x } t)\ 2 dx. (4.1.11) 

The role of this quantity, which should be positive by definition, is fundamental in the 
Zakharov's theory of the stability [[T7], 78]. From (|4.1.1| ) and (|4.1.11|) , we get 



d 



(\x\ 2 ) = - f((f>*x- V(j)-(j)x- V(jf)dx 

= 2 Im J V0dx, (4.1.12) 



where the superscript * means the complex conjugate. By differentiating Q4.1.12] ) with 
respect to t and using Eq. ( |4.1.1|) again, we obtain [M 



dt 2 



\x\ 2 ) =2H- j (2V{x) + x ■ VV(x))\(f)\ 2 dx - ^(D - 2) ||0||t (4.1.13) 



Substitution of (|4.1.2|) into ( |4.1.13|) yields the equation of motion for (|a;| 2 ), 



dt 

where 



f-(\ x \ 2 ) = 2H-AA 2 (\x\ 2 )-f(t), (4.1.14) 



/(f) = 4 J2(u 2 n - A 2 )\\x n m\l + \{D - 2) U(t)\\i (4-1.15) 

n=l Z 
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with A being denned by ( |4.1.8| ). Note that f(t) > 0, if D > 2. A general solution of the 
differential equation (|4.1.14| ) is 



f(u) sin[2A(t - u)]du, 



where A and 9q are constants. Without loss of generality, we assume that 

A > 0, < 6 < 2tt. 

For the later convenience, we define a function, 1(9), 

H 



1(9) 



+ A sin 9, 



with 



Then, ( |4.1.16| ) becomes 



2A 2 
9 = 2At + 9 . 

1 /"* 



A 2 ) = IW)) -—J f(u) sin[2A(t - u)]du. 



(4.1.16) 
(4.1.17) 

(4.1.18) 
(4.1.19) 

(4.1.20) 



The constants A and 9q are related to the initial conditions of the wavefunction. Using 
( pTTTTTD and (|4.1.16|) in ( ^PO]) , we get 



F t=o 



l(9 



x\ 2 \4>o(x)\ 2 dx 
H 



2A 2 



+ A sin 9 , 



where 



4>o(x) = <p(x,t = 0). 
Differentiation of ( [4.1.16 ) with respect to t gives 

^-(|x| 2 ) = 2AAcos(2At + # ) - /*/(«) cos[2A(t-u)]du. 
dt jo 

Then, from (|4.1.12|) and ( [4.1. 23|) , we get 



(4.1.21) 
(4.1.22) 

(4.1.23) 



dt 



(kl 2 ) 



t=0 



t y (00 s • Wo -(f>o%- V0q) da; 
= 2Av4cos6> . 
Solving ( [1.1.211 ) and ( |4.1.24| ), we obtain 



(4.1.24) 



A 



(<M 2 >l--^)*+(^>l 2 > 



1 d 



t=0/ 



1/2 



(4.1.25) 
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sin #0 
cos 6> 




(4.1.26) 
(4.1.27) 



In the following, we analyze the solution (|4.1.16|) or equivalently ( fi.1.20 ) in the case 
that the spatial dimension is larger than or equal to two, namely D > 2. 
1) H < case 

We consider the case that the total energy of the system H is less than or equal to zero. 
In this case, we show that the wavefunction surely collapses in a finite time. 

From (|4.1.21|) , l(9 ) should be positive because of the positivity of (\x\ 2 ) defined by 
( |4.1.11|) , and therefore the constant A should satisfy 



A > 



H 



2A 2 



H 



2A 2 



(4.1.28) 



When ([4.1.28|) is satisfied, there exist two zero points of 1(9) which lie in < 9 < ix. We 
write them as 0+ and 9- (0 < 9- < n/2 < 9 + < ir). Explicitly, 9± are given by 



+ 
9_ 



n — Arcsin ( 
-H 



, -) 

\2A 2 AJ 



Arcsin 



2A 2 A 



(4.1.29) 
(4.1.30) 



Here and hereafter, the region (branch) of Arcsinx is defined to be [— tt/2, +7r/2]. Then, 
the constant 9q should be 

9_ < 9 < 6+. (4.1.31) 

Since the function f(t) defined by ( |4.1.15| ) is larger than or equal to zero when D > 2, 
and sin[2A(t — u)] appearing in ([4 .1.1 6|) is also larger than or equal to zero for < t < 
we obtain an inequality, 

(\x\ 2 )<l(9(t)), (4.1.32) 

for < t < £r or equivalently 9q < 9 < 9q + n. The function 1(9 '(£)) becomes surely 
negative for some value of t which lies in < t < £r, when ([4. 1 . 28| ) and ( [4.1 . 31|) are 
satisfied. Then, from the inequality (|4.1.32|), (|a;| 2 ) should also be negative for some value 



of t, < t < This contradicts to the definition of (\x\ 2 ) ( |4.1.11| ) and implies the 
development of the singularity, characterized by the divergence of || V0||2, in a finite time 
defined as t . In fact, from the Cauchy-Schwarz inequality ( |4.1.9| ), we see that || V0||2 tends 
to infinity when (|x| 2 ) goes to zero. It is known that the instant to ma y differ from the 
moment at which (|x| 2 ) vanishes, defined as t c , and that to is l ess than or equal to t c , in 



general [pq| . We call t c the collapse-time in this paper. 

We can show that sup^gpo \4>(x, t) \ tends to infinity when || V0||2 diverges. This problem 
has been discussed also in Refs. |79| , ^5| . For g — — 1, the energy ( [4.1.4|) with the harmonic 
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potential ( |4.1.2| ) is written as 

H = /(W + £ 

\ n=l 



«H a - -101*1 da 



D 



l|V0|| 2 + £* 



'112 



n=l 



Then, the following inequalities hold, 



D 



-2iJ + 2||V0||2 + 2£ 



v n\\ x n<P\\2 



n=l 

> -2H + 2\\Vm + 2A 2 (\x\ 2 ) 

> — 2^T + 2|| V0||1 + 2(AZ?A^/2) 2 ||V0||; 



(4.1.33) 



(4.1.34) 



In deriving the third line of ( 4.1.34 ), we have used ( |4.1.9|) . Further, combining an inequality 



\< sup \</>{x,t)\) U\\l= sup \<f>(x,t)\) N, 



(4.1.35) 



with ( |4.1.34| ), we find that sup^pjo \(p(x,t)\ goes to infinity when ||V0||2 diverges. 

From the inequality ( |4.1.32| ), we see that (|x| 2 ) vanishes before l(9(t)) reaches zero at 
9 = 9 + . Thus, by using ( |4.1.19|) , we can obtain an upper bound of the collapse-time, t CjM , 



When 9 < n/2, namely ^(\x\ 2 ) t _ n > 0, from (|4.1.26| ) we get 



t=o 



9q = Arcsin 



' (\x\ 2 )\t- 
A 



— V 

2A 2 A 



Substituting ( [4.1. 29|) and (|4.1.37| ) into ([4.1.36|) , we have 



t-c,M 



1 

2A 
1 

2A 



7r — Arcsin 



\2A 2 AJ 



Arcsin 



\x\ 2 )\ t =o H \ 
2A 2 A 



A 



where £ is defined as 

c = 



tt/2 + Arcsin £), 



2A dt 



t=0/ 



X 



2A dt 



\A 2 ) 



-A 



(N 2 )Uo (|x| 2 )| t=0 



A 2 



J_ d_ 

2A d7 



<M 2 > 



(4.1.36) 



(4.1.37) 



(4.1.38) 



t=tv 



+ ^((N 2 )U-^) 



(4.1.39) 
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In deriving Q4.1.3S 



), we have used a formula, 



Arcsin 77 — Arcsin ( 



Arcsin [r]Jl - ( 2 - (Jl - r} 2 , (4.1.40) 



for < 77 < 1 and < ( < 1- 



When 9 > vr/2 namely 4:{\x\ 2 ) 



< 0, 6*o becomes 
t=o ~ 



instead of (|4.1.37| ). In spite of this change, we obtain t c M just the same as (|4.1.38|) . 



2) H > case 

For the H < case, we have shown that the wavefunction surely collapses in a finite 
time, regardless of the initial condition on the wavefunction 4>o(x). Here we prove that 
even when the total energy H is positive, the collapse of wavefunction occurs in a finite 
time for a certain class of the initial conditions. We investigate such conditions according 
to the value of 9 , as follows. 

We first consider the case that < 9q < tt/2, which is equivalent to 



> 0. (4.1.42) 



t=o 



As mentioned in the H < case, f(t) defined by ( (4.1.15[ ) is larger than or equal to zero 



when D > 2, and sin[2A(t — u)} appearing in (|4.1.16|) is also larger than or equal to zero 
for < t < or equivalently 9 < 9 < 9 + tt. Then, if 

l(9 + n)<0, (4.1.43) 

1(9 (t)) becomes negative for some value of t in an interval < t < and thus due to 
( [4.1 .32| ) the collapse of the wavefunction happens. By use of ( [4.1 .21| ) and 



l(9 + 7r) = H/A 2 -l(9 ), (4.1.44) 
the condition (|4.1.43|) is rewritten as 

(|x| 2 )|* = o>i//A 2 . (4.1.45) 



When the condition (|4.1.45| ) is satisfied, there exist zero points of 1(9), of which we denote 



the positive and minimum one by 9*. It is easily confirmed that 9* lies in n < 9* < 3tt/2. 
Considering < 9 < 7r/2 and using (|4.1.26|) , we get 



A 



9 = Arcsin ( - — 1 . (4.1.46) 
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Similarly, since 1(9*) = and n < 9* < 3tt/2, we have 

9* = 7r + Arcsin 



(—) 



(4.1.47) 



Then, from (|4.1.46[ ) and ( [4.1.47| ), we obtain the upper bound of the collapse-time £ Ci m, 



<-c,M 



= ^ (tt/2 + Arcsin 0, 

with £ being defined by ( |4.1.39| ). 

Next, we consider the case that tt/2 < 9 < 3ir/2, namely 



(4.1.48) 



< 0. 



In this case, if 



A > 



H 

2A 2 ' 



(4.1.49) 



(4.1.50) 



l(9(t)) becomes negative for some value of t which lies in < t < Then, if ( [4.1. 50|) is 
satisfied, we see from ( |4.1.32| ) that collapses in a finite time, and the upper bound of the 



collapse time t c> M is calculated in the same way as ([4.1.48 



We further investigate the collapse-condition ( [4.1.501) . Substituting (|4.1.25| ) into (|4.1.50j ), 
we have 



2A dt 



(kl 2 ) 



> (\x\ 



\t=0 



t=0/ 



H 

A2 



<N 2 )I 



i=0 



(4.1.51) 



When if/A 2 < (|x| 2 )| t= o which is equivalent to the collapse-condition for the < 9 < n/2 
case ( [4.1.45| ), the right hand side of ( |4.1.51| ) is negative. Then, regardless of the value of 



i ^(|x| 2 ) (|¥X5TD holds. When H/A 2 > (|x| 2 )|i =0 , the right hand side of (gXH]) is 



positive. By taking the square roots of both sides of (|4.1.51|) and noting that ^(|a;| 2 ) 
is negative, we get 



t=o 



< 



-2^(\x\ 2 )\ t=0 (H-A 2 (\x\ 2 )\ t=0 ). 



(4.1.52) 



t=o 



The condition (|4.1.52 ) can be regarded as an extention of the one obtained by Weinstein 
in the case (iii) of Theorem 4.2 in Ref. [BH], where the "free" or the conventional (i.e. 
without external potential) nonlinear Schrodinger field is considered. In fact, if v n = for 
all n, ( 4.1.52| ) becomes the same (apart from the discord of the factor, 2) as the one for 
the free nonlinear Schrodinger (NLS) field [80]. We note that, if 

v n — v (for all n), (4.1.53) 
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we have 

H - A 2 (|x| 2 )|, =0 = || VMl ~ Ml/2, (4.1.54) 
and thus the collapse-condition ( [4.1.52|) becomes formally the same as in the free NLS field 



case. 

Finally, we consider the case that 37r/2 < 9 < 2tt, namely 



> 0. (4.1.55) 



From ( |4.1.21| ), l(9 ) should be positive, yielding that the function 1(9) in ( |4.1.18| ) is always 
positive ior 9 < 9 < 9 Q + n or equivalently < t < Then, we cannot see from ( [4.1. 32|) 



whether (|x| 2 ) becomes negative for < t < Further, for t > the inequality ( [4.1.32| ) 
does not hold any more, because the integral 

f f(u) sin[2A(t - «)]d«, (4.1.56) 
Jo 

appearing in the right hand side of (|4.1.16| ) is not necessarily positive for t > Accord- 
ingly, we cannot see whether (|x| 2 ) becomes negative for t > £r, either. Therefore, in the 
case that 37r/2 < 9q < 2tt, we cannot conclude that the wavefunction collapses in a finite 
time. 

It is interesting to consider a special case that the space is two-dimensional (D = 2), 
and the constants {u„} satisfy ([4.1.53|) . In this case, f(t) defined by Q4.1.15|) is identically 
zero, and ( |4.1.16j ) reduces to 

H 

2^2 



(|x| 2 ) = l(6(t)) = — + Asm(2vt + 9 ). (4.1.57) 



Then, we can trace the exact time-evolution of (|x| 2 ). If the collapse of happens which 
can be proved in the same manner as in the previous cases, the collapse-time t c is equal to 
its upper bound t CjM . 

To conclude this sub-section, some comments are in order. First, what we have pre- 
sented in this section is a sufficient condition for the collapse. There remains a possibility 
that, even when the collapse-condition in this section is not satisfied, the expectation value 
(\x\ 2 ) may go to zero after several damped oscillations, which is due to the effect of the neg- 
ative non-homogeneous term —f(t) (D > 2) in the ordinary differential equation ( |4. 1 ■ 14|) . 



In this case, the upper bound of the collapse-time t c ^ M ( [4.1.48| ) will be replaced by larger 
ones. A quantitative analysis of f(t) is required to make clear this damped-oscillating 
phenomenon, which will be a future problem. 

Second, different from the discussion in the repulsive case (g = +1), we have not 
considered the case D = 1 here. If D — 1, the positivity of f(t) defined by ( [4.1. 15| ) is not 
guaranteed, and thus we cannot examine the stability of the wavefunction by using the 
above mentioned method. The one-dimensional and self-focusing NLS equation with no 
external potential is integrable and has stable solition solutions []73]. The external harmonic 



potentials destroy the integrability. However, since the equation has no singularity at the 
origin, we conjecture that the effect of the potentials is only a deformation of solitary 
waves. 
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4.2 Collapse of the condensate 



We predict the collapse of the Bose-Einstein condensate in a magnetic trap where the 
atoms have a negative s-wave scattering length [41]— 1|43|. We prove that the singularity of 
wavefunction emerges in a finite time even when the total energy of the system is positive. 
In addition, we present a refined formula for a critical number of atoms above which the 
collapse of the condensate occurs. This number for 7 Li atoms can be the same as the one 
in the recent experiment. 



4.2.1 Time- evolution of the condensate 

We consider the Bose-Einstein condensate confined by magnetic fields and investigate the 
time evolution of wavefunction when the interactions between atoms are effectively attrac- 
tive. Since the magnetic trap has an axial symmetry and the system is near the ground 
state, we assume that the wavefunction is axially symmetric. The axis of the symmetry 
is chosen to be the z-axis. Let t denote time and r± the radius of the projection of the 
position vector r on the xy-plane. Time-evolution of macroscopic wavefunction ^f(r±, z,t) 
of the condensate at very low temperature is described by the Gross-Pitaevskii equation 
for the system |7D[-|72|, 



. t 3* n 2 l d ( h 2 d 2 ^> l 22t i 2 2t 47rftV Tl2T , Ani . 

Here m is the atomic mass, u z and u± are the trap (angular) frequencies along the z-axis 
and in the xy-jAane, respectively, and a is the s-wave scattering length. We take a to be 
negative. By introducing a characteristic length ro, 

( h \ 1/2 

we prepare dimensionless variables, 

P = r±/r , £ = z/r , r = u±t, ip = r^ 2 ^. (4.2.3) 
In terms of these variables, Eq. (|4.2.1| ) is written in the dimensionless form 

where 

5 = Uz/uj^, c = 87ia/r . (4.2.5) 

Note that the constant c is dimensionless and negative. In the terminology of the soliton 
theory, we refer to Eq. fl4.2.4Q as axially symmetric attractive nonlinear Schrodinger (NLS) 
equation with harmonic potentials. The NLS equation ( 4.2.4 ) has two integrals of motion, 



OO POO 



h = 2n dU pdp\W, (4-2.6) 



o 
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/oo rco 
d£ / pdp 
-oo Jo 



(4.2.7) 



The total number of atoms N and the total energy E are related to the constants Ii and 
h as 

N = h, E = hu ± I 2 . (4.2.8) 



By definition, I\ is positive. We define 

/oo poo 
ad P dp(p 2 + eM 
-oo JO 



(4.2.9) 



which measures the geometrical extent of the wavefunction ip and plays a key role in the 
Zakharov's theory [[T?], |78| . 

Differentiating (|4.2.9|) by r twice and using ( |4.2.4| ), we have 



dr 2 



(v 2 ) 



411V) -/(t), 



where Q and /(r) stand for 



il = min(l, 5), 

/oo poo 
dd pd P [2(5 2 -i)em 2 + \c 
-co JO 

/co rco 
dd pdp[2(l-5Vh/f + 
-oo Jo 

A differential equation ( [4.2. 10|) for (rj 2 ) is readily solved to give 



4 ] (for 5 > 1), 
4 ] (for < 5 < 1). 



(rT°> =Asin(2ftT + # ) + — / 2 



1 

2Q Jo 



f(u) sin[2fi(r — u)]du. 



(4.2.10) 

(4.2.11) 
(4.2.12) 

(4.2.13) 
(4.2.14) 



Constants A and 8q are to be determined by the initial conditions on ip. Without loss of 
generality, we assume A to be positive. 

We first consider the case that I 2 and therefore the total energy of the condensate E 
are negative. Because (rj 2 ) is a positive quantity, the right hand side of ( |4. 2. 14| ) should 
also be positive for r = 0. Accordingly, 9 in (f4.2.14|) must be < 6> < ix. Let us define a 
function, G(r), as 

G(t) = A sin(2fir + 6 ) + -| J 2 . (4.2. 15) 

Since the function /(r) is always positive and sin[2r2(r — u)] in ( f4.2.14| ) is positive for 
< t < ir/2tt, the last term in (|4.2.14j ) is negative for < r < tt/2Q. We thus obtain an 
inequality, 

(V 2 ) < G(t), (4.2.16) 

for < r < tt/2Q. From the definition ( |4.2.15| ), we see that G(r) becomes negative for 
some value of r which lies in < r < ir/2Q, if < 6q < n. Then, from the inequality 
( [4.2. 16| ), (rj 2 ) should become negative for some value of r, < r < tt/2Q. This contradicts 
to the definition of (if) ( [4.2.9|) , which implies the development of the singularity. Thus, 
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when the total energy of the condensate E is negative, the wavefunction ip surely collapses 
in a finite time, regardless of the initial condition on ip. 

Next, we prove that, even when the total energy E is positive, the collapse of wavefunc- 
tion occurs in a finite time for a certain class of the initial conditions on ip. Actually, in 
the experiments of the Bose-Einstein condensation under a magnetic trap, the total energy 
E and therefore I2 are positive. We assume the initial shape of %fj to be gaussian, 



/ 2^1/27- \ 1/2 

*{p,lt = o) = [q ex P [-( g y + sm/ii 



(4.2.17) 



where q and k are positive parameters. From Eqs. (|4.2.4j) , ( |4.2.9 ) and ( |4.2. 1 7| ) , we can 
show that 

^-(rf) = (forr = 0), (4.2.18) 
dr 

which yields 9 = n/2 or 9 Q = 3n/2 in ( |4.2.14j ). Correspondingly, the solution ( |4.2. 14|) now 
yields 



( V 2 ) = ( V 2 ) + = Acos(2fir) + Jq 



f(u) sin[2f2(r — u)]du, 



(4.2.19) 



when 9 = 7r/2, and 

( v 2 ) = ( V 2 )_ = -Acos(2Qr) + A/ 2 _ _L ^ T /(m ) sin [ 2 fi(r - it)] du, (4.2.20) 

when 6*o = 37r/2. Substituting the initial shape ( |4.2.17| ) into (rj 2 )\ T=0 and I2, we get 

(4.2.21) 



(r ) 2 )\ T=0 = (2q- 2 + 6- 1 k- 2 )I ll 



\{q 2 + q - 2 ) + - A {k 2 + k- 2 ) 
From Eq. ([4.2.22 ), we see that I2 is positive when 

h < h >a (q, k). 

Here 



/l + l6^72- Jl - 



A,a(g, k) = 



47T 3/2 



\c\q 2 k8 1 / 2 

On the other hand, Eqs. ( |4.2.21| ) and ( |4.2.22| ) give 



[2(q 2 + q- 2 ) + 5(k 2 + k- 2 )]. 



(V 2 )\r=0 ~ 



\c\q 2 k5 1 / 2 

87r 3/2 2 h(h-h,o(q,k)), 



where h,b(q, k) is defined as 



47r 3/2 

J 1>H9' k ) = l-l-2 fc gl/2 ^(g' fc )' 



|c|g 



(4.2.22) 

(4.2.23) 
(4.2.24) 

(4.2.25) 
(4.2.26) 
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with 

g{q, k) = 2[q 2 + (1 - 2tt 2 )q~ 2 ] + 5[k 2 + (1 - 2n 2 /5 2 )k~ 2 }. (4.2.27) 

Depending on whether I\ is larger or smaller than h^iq, k), we can determine which of 
(r] 2 )± should be chosen as the solution of Eq. ( [4.2.101 ). When Ii > I±fi{q, k) (Ii < h,b(q, k)), 
{rf) + ({if)-) is the solution of Eq. ( |4.2.10| ) . Remark that the function g(q, k) may take both 
positive and negative values, because 1 — 2Q 2 = — 1 < for 5 > 1 and 1 — 2(Q/5) 2 = — 1 < 
for < 5 < 1. 

In the case that g(q,k) > and therefore Ii t b(q,k) > 0, if I\ < Ii t b(q,k), (i] 2 )- is an 



appropriate solution of Eq. ( |4.2.10| ) and satisfies 



(v 2 )- 



< 



-Acos(2fir) + — 1 2 , 



(4.2.28) 



for < r < 7r/2fi. The right hand side of the inequality (f4.2.28| ) is always positive because 



.4 



1 T 



{rj 2 )\ T= o < 2I 2 /fl 2 . Then, we cannot see from (|4.2.28| ) whether (r/ 2 )- surely 
becomes negative for < r < tt/2Q. Further, for r > ir/2Q, the inequality ( [4.2.28| ) does 
not hold any more, because the integral, 



f(u) sin[2f2(r — u)]du, 



(4.2.29) 



is not necessarily positive for r > tt/2Q. Accordingly, we cannot see whether (?7 2 )_ becomes 
negative for r > tt/2Q, either. Therefore, when 1\ < Ii^{q,k), we do not conclude that 
the wavefunction ip collapses in a finite time. This conclusion is consistent with that in 
the case 3n/2 < 9 < 2tt in the sub-section [4.1. 2| . On the other hand, if Ji > Ii t b(q,k), 
the behavior of (r] 2 ) may change drastically. As the solution of (|4.2.10|) , (i] 2 ) + defined by 
( |4.2. 19[ ) is appropriate and satisfies 



( V 2 ) + < Acos(2Qt) + —I 2 , 



for < r < vr/2fi. Using (|4.2.19|) and (|4.2.25|) , we have 



A=(V 2 ) \ T =»-- 2 h = -— 2 



g(q,k)h + 



\c\q 



Suppose that 



a > 2i 2 /n 2 



2 k^l 2 



87T 3 / 2 fi 2 



(4.2.30) 



(4.2.31) 



(4.2.32) 



which is possible as seen from ( f4.2.31| ). When the inequality ( |4.2.32j ) holds, the right hand 
side of ( |4.2.30| ) becomes negative for some value of r which lies < r < 7r/2f2. Then, 
from ([4.2.30|) , (t] 2 )+ also becomes negative for some value of r, < r < tt/2Q. This is 
contradictory to the positivity of (i] 2 )+, implying that the singularity of the wavefunction 
emerges in a finite time. In other words, the wavefunction ip collapses in a finite time. Let 
us rewrite the condition ( |4.2.32| ) so as to make its physical significance clear. Substituting 
( ^2~22l) and (|4.2.31|) into ( gX3|) , we get 



h > h, c {q,k), 



(4.2.33) 
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where 

4^3/2 

h ^ k) ^Wk^ Kq ' kl (4 ' 2 ' 34) 

with 

%, k) = 2[q 2 + (1 - tt 2 )q~ 2 } + 6[k 2 + (1 - n 2 /5 2 )kr 2 }. (4.2.35) 

Contrary to g(q,k), the function h(q,k) is always positive. From ( |4.2.24j ), ( f4.2.26| ) and 
Q4.2-34|) , we find that 



h,b < h,c < h,a (for all q and k). (4.2.36) 

Therefore, if the initial shape of ip is gaussian (|4.2.17|) , the collapse occurs when I\ > 
Ji jC (g, k). The inequality ( |4.2.36| ) assures the consistency of the analysis. 



In the case that g(q,k) < 0, and therefore Ii t b(q,k) < 0, {rf) + (|4 .2.1 9|) is the only 
solution of ( |4.2.10| ) because I\ should be positive. Except for this, the discussion for the 



collapse goes in the same manner for the h^(l, k) > case. Thus, we conclude that the 
solution of the NLS equation ( f4.2.4| ) with the initial condition ( f4.2.17| ) surely collapses in 



a finite time for I\ > Ii jC (q, k) even when the total energy is positive. 

The above result predicts an interesting phenomenon, the collapse of the Bose-Einstein 
condensate, for the assembly of bosonic atoms with a negative s-wave scattering length. 
Using (|4.2.5|) , Q4.2.8|) and ( |4.2.34|) in (|4.2.33| ), we obtain a critical number of atoms N c , 



N c (q,k)=N ^^, (4.2.37) 

with 

N = (7r/5) 1 / 2 r /(2|a|). (4.2.38) 

The collapse of the wavefunction occurs when the number of the trapped atoms, N, exceeds 
N c (q, k). The formula fl4.2.37|) gives the critical number of atoms iV c as a function of two 
parameters q and k in ( [4.2. 17| ). The extents of the initial wavefunction are proportional to 
q^ 1 for the p-direction and k' 1 for the ^-direction. When q and k equal 1, the initial wave- 
function ( |4 . 2 . 1 7| ) represents the ground state for the Schrodinger equation under harmonic 
potentials, 

It is interesting that the formula includes another anisotropy effect of the magnetic trap 
through Q = min(l, 5 = uj z /ujj_). 

4.2.2 Application to 7 Li system 

We estimate the critical number for the assembly of 7 Li atoms which have a negative s-wave 
scattering length. In order to apply the above discussion to this system, we first improve 
quantitatively the evaluation for (r] 2 ) + ( |4.2.30|) as follows. If there exists a constant, /, 
such that < / < /(r), the inequality ( |4.2.30| ) can be replaced by 

<" 2> + £ { A + w?) cos(2f!T) + ii h ~ w?) • (42 ' 40) 
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for < t < n/2Q. Then, the condition for the collapse of ip is modified as 



2 f 



fl 2 z 2ft 2 ' 



(4.2.41) 



instead of ( |4.2.32| ). Here, we assume that the collapse of the wavefunction occurs without 
oscillations, which means that (t] 2 ) goes to zero monotonously and the integral, 



/oo roo 
-oo JO 



(4.2.42) 



is a monotone increasing function of a (scaled) time r. Under this assumption, from ( [4.2. 12| ) 
or ( [4.2.13|) , we can estimate / as 



/oo roo 
d£ / pdp\ip(p,£,T = 
-oo JO 

Using the initial condition on ip ( [4 . 2 . 1 7| ) in ( [4.2. 43[ ), we get 

f = I 2 kq 2 /N . 



(4.2.43) 



(4.2.44) 



Thus, substituting ( |4.2.22| ), ( |4.2.31| ) and ( |4.2.44[ ) into ( |4.2.41| ), we arrive at a renewed 
critical number of atoms, iV C)new , 



N, 



c,new -^0 



^new(O') ^) 

q 2 k 



where 



Kcw(q,k) = -h(q, k). 



(4.2.45) 



(4.2.46) 



Apparently, iV Cjnew is smaller than N c . 

In the experiment for 7 Li | 12]| , the values of frequencies uj z and were uj z /2tt = 117 Hz 
and co>_|_/27r = 163 Hz, which gives 5 = 0.718. The s-wave scattering length of 7 Li was 
observed to be a = — 27.3a (ao : Bohr radius) P§| . Substituting these experimental 
values into (|4.2.38j) and (|4.2.46|) , we get 



iVc,new(g, k)/N = 4AT 1 (1 + 0.485g" 4 )/3 + 0A79q- 2 k, 



(4.2.47) 



with iVo = 1520. 

To fix a relation between q and k, we use the minimal condition for I2 ( [4.2.22| ). From 
dgh = and dkh = 0, we have 



q 2 - q - 2 = S(k 2 - k 



-2\ 



which yields q as a function of k, 



q = Q{k) 



5{k 2 - k- 2 )/2 + ^l + 5 2 (k 2 -k- 2 ) 2 /A 



1/2 



(4.2.48) 



(4.2.49) 
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Figure 4.2.1: N c ^ new /N as a function of k for the case that 5 = 0.718 [0. 



Substitution of 5 = 0.718 into ( [4.2.490 gives 



q = Q(k) = 



0.359(V - fc ) + t/1 + 0.129(A; 2 - jfe 



-2^2 



,1/2 



(4.2.50) 



We plot N„{Q(k) } k)/N Q as a function of fc for 5 = 0.718 (Fig. gjj). The refined 
critical number N c ^ new (Q(k), k) is a monotone decreasing function of fc. This behavior is 
reasonable: For smaller (larger) k, the wavefunction is more spread (confined), and such 
configuration imposes more (less) severe restriction on the collapse. In what follows, we 
consider the case k > 1 since the gaussian wavefunction ( [4.2. 17|) with q = k = 1 does not 
include the effect of attractive inter-atomic interaction which makes the exact ground state 
more confined. This has been considered in the sub-section 3.1.3 . For instance, in the case 
of almost spherical trap, the result ( ^.1.66| ) gives k = 1.59 and q = 1.45. 

If we set k = 1.20, Q(k) and N rnpvir (Q(k),k) are calculated as Q(1.20) = 1.14 and 



iVc,new(<3(1.20),1.20) = 2840. In the case that k = 2.20, Q(k) and N c>nevr (Q(k),k) are 
estimated as Q = 1.90 and iV C)new = 1400. The latter value of iV c>new agrees well with the 
one in the experiment, which is about 10 3 [13|. It is also consistent with the theoretical 
ones obtained by different approaches, which are (1 ~ 3) x 10 3 pi]- ||40|| , [^5|-[pT[ . 



4.3 Summary 

We have studied the dynamical properties of Bose-Einstein condensates. In Sec. W7L 



as 



a preliminary to the analysis of the collapse of the condensate, we have considered the 
stability of the wavefunction of the D-dimensional nonlinear Schrodinger equation ( |4.1.1[ ) 
with harmonic potential terms ( |4.1.2|) , for both repulsive (g = +1) and attractive (g = — 1) 
cases. This is an extention of the Zakharov's theory (77], Q . In the repulsive case, we have 
shown that the solution of Eq. (|4.1.1|) is absolutely stable for any spatial dimension when 
confined in the harmonic potential. In the attractive case for D > 2, by solving the equation 
of motion for the expectation value of the square of the radius ( |4.1.14j ), we have investigated 



the time-evolution of the wavefunction (p. When the total energy of the system H is equal 
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to or less than zero, the wavefunction <p surely collapses in a finite time, regardless of the 
initial condition. On the other hand, when H > 0, we have obtained conditions that the 
collapse of the wavefunction occurs, according to the value of the constant 9 appearing 
in ( 4.1.16 ). When < 6q < n/2 which is equivalent to ( 4.1.42 ), the collapse-condition is 
obtained as ( |4.1.45|) . When ti/2 < 9 < 37r/2, equivalent to (|4.1.49| ), the collapse-condition 
is ( |4.1.52|) . This result is an extention of Weinstein's where the nonlinear Schrodinger 
(NLS) field without external potential, is considered ||30| . We have also obtained an upper 
bound of the collapse-time of the wavefunction £ Ci m, represented by the formula ( |4.1.38| ) 
with ( 4.1.39| ). In a special case that D = 2 and v n = u, since the function f(t) defined by 
( |4. 1 . 15[ ) is identically zero, the time-evolution of (\x\ 2 ) can be traced exactly by ( (4.1.57[ ), 
and therefore the collapse-time t c is just equal to its upper bound t CjM . Improvements of 
the collapse-conditions and the upper bound of the collapse-time require a more detailed 
analysis of f{t). We, however, emphasize that the collapse time is the order of \jv. 

In Sec. [4.2| , we have considered the instability of the Bose-Einstein condensate with ef- 
fectively attractive inter-atomic interactions under a magnetic trap. The Gross-Pitaevskii 
equation (|4.2.1|) for the system is essentially the same as the attractive nonlinear Schrodinger 
equation with harmonic potential terms ( |4.2.4| ). First, we have shown that the collapse of 
the wavefunction surely happens in a finite time when the total energy of the condensate E 
is negative. Second, by assuming the initial condition on the wavefunction to be gaussian 
( [4.2. 17|) with two parameters q and k, we have proved that the collapse can happen even 
when E is positive. This situation corresponds to the experiments on the Bose-Einstein 
condensate. We have also proved that the condensate collapses when the number of atoms 
N exceeds N c ( |4.2.37|) . Further, we have presented a refined formula for the critical num- 
ber iV C)new ( |4.2.45| ), which is smaller than N c . Using a relation between q and k which is 
obtained from the minimal condition for the energy (|4.2.49| ), we have plotted N c ^ new /N 
as a function of k for 5 = 0.718 (Fig. 4.2.1 ). Our estimated number for the assembly 
of 7 Li, which is about 1400 for q = 1.90 and k = 2.20, agrees well with the one in the 
experiment [[13]. We have simplified the analysis by assuming (i) axial symmetry of the 
system, (ii) zero temperature and (iii) s-wave scattering. All the assumptions are more or 
less consistent in the sense that we only consider the properties of the dilute condensate at 
the very low energy region. In addition, we have not taken account (iv) three-body colli- 
sions and (v) the short-range repulsive part of the inter-atomic potential, which may play 
important roles at the final state of the collapse. Detail examinations of those neglected 
effects are left for future studies. On the other hand, it is extremely interesting to observe 
the time-development of the wavefunction near the collapse. The asymptotic analysis has 
been given in 
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Chapter 5 

Summary and Concluding Remarks 



We have reviewed recent theoretical research on the properties of Bose-Einstein condensa- 
tion, in particular, static and dynamical stabilities of Bose-Einstein condensates confined 
in magnetic trap. It has been shown that the variational approach and the analysis of the 
time-dependent nonlinear Schrodinger equation give significant information on the stability 
of the condensates. 

Bose-Einstein condensation was discovered in liquid helium about 70 years ago, and 
has been "rediscovered" in alkali atomic gases confined in trap, due to the developments 
of technology in atomic physics. There are many advantages in experiments using atomic 
vapors under magnetic trap, which are summarized as follows. 

1) The atom-atom interactions are weak (the s-wave scattering length is about 10~ 8 m, 
whereas at the required densities the inter-atomic spacing is about 10 _6 m, for the 87 Rb 
case Therefore, the purely quantum-statistical nature of the BEC transition can be 
captured. 

2) We can choose various species of bosonic (or fermionic) atoms and prepare even their 
mixtures with different number of each atoms. 

3) We have controllable parameters such as density, temperature and the frequencies of 
the magnetic trap. In particular, the number of atoms and the temperature are now 
independent. This offers an idealistic situation to study the quantum statistical mechanics. 

4) The interaction strength can be varied by use of the Feshbach resonance H88 , S9| . The 



s-wave scattering length depends on the applied magnetic field and even the sign of it 
changes. 

These new aspects give rich dynamical behaviors of the condensate at microscopic, meso- 
scopic and macroscopic levels. Understanding of those systems is fundamental in quan- 
tum statistical mechanics, condensed matter physics and atom optics. The trapped Bose- 
Einstein condensates may serve as a source of the coherent matter waves. In atom optics, 
we expect many applications such as atom holography P0| , pi], atom interferometer |J2 



and atom laser [P5|-|P3|. In principle, Bose-Einstein condensation is expected to occur for 



all bosonic atoms and molecules with and without charges, and thus we believe that those 
studies will continue to be interesting in coming years. 
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Appendix A 

Bose-Einstein Condensation of a Free 
Boson Gas Confined in Harmonic 
Potentials 



We discuss the Bose-Einstein condensation of a free gas under harmonic potentials. The 
iV-body Hamiltonian (first quantization) is 



H 



N 
N 



E ( 2 i 2 i 2 \ i ' ' l ( 2 2 i 22, 



3=1 



m 
2~ 



2 2 



(A.I) 



where = (Pj X ,Pj y ,Pjz) is the momentum operator of particle-j, = —ikd/drj. Since 
there is no interaction between particles, the eigenstates are expressed as 



N 



tf(ri,r 2j - ■ -,r N ) = H^rf, 

3=1 



(A.2) 



and a single-particle state <p(r) satisfies 



h 2 ( d 2 



d 2 d 2 

+ T^T + 



2m y<9x 2 <9y 2 <9z 2 
The eigenvalues of ( |A.3| ) is well-known; 



e0 r* 



(A.3) 



-n x n y n z 



(A.4) 



where n x , n y , n z = 0, 1, 2, 
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We adopt the grand canonical ensemble. The total particle number N and the total 
energy E are given by 



N = Yl [ eX P f 3 { e n x n y n z ~ I*) ~ 1 

Tlx i^wi^z 



-i 



E 



E 



exp 



ft {^n x n y n z l^j 1 



(A.5) 
(A.6) 



where (3 = (k-sT)~ l and \i is the chemical potential. For later discussions, we shift the 
chemical potential as /i — h(u x + uj v + oj z ) /2 — > //. At very low temperature, ( |A.5| ) is 
written in the following form, 

N = N Q + N 1 , (A.7) 



where 



N 



l/(e- 



1 



[exp (5 (n x huj x + riyhujy + n z hu z ) — 1] 1 . 

n x ,n y ,n z =£0 



(A.8) 
(A.9) 



Discussions, here and in what follows, are essentially the same as those for a free boson gas 
in a box. The Bose-Einstein condensation is the situation where iVo becomes macroscopic, 
that is, N ~ O(N). Equation flOj) implies n < and -fin ~ 0(1/N). We have set 
H = in ( |A.9|) . Therefore, Ni in ( |A.9|) gives the maximum of the contribution from the 
excited states at low temperature. 

We replace the summations in ( |A.9| ) by an integral over the single-particle energy e 
with the density of states D(e), 



D(e)de 



-de. 



(A.10) 



The formula (|A.10|) can be derived as follows. We estimate the number of the states, N(e), 
whose energies are smaller than e, which is equivalent to the number of positive integer sets 
{n x ,n y ,n z } satisfying < n x hu x + n y hu y + n z hu z < e. Geometrically, this corresponds 
to a volume, (e/huj x )(e/hLjy)(e/hLJ z )/6. Therefore, we obtain N(e) = ^e 3 /(hu) 3 , where 
Cu 3 = uj x uj y uj z . A relation D(e) = dN(e)/de gives (|A.10|) . 

Using ( A.10|) to rewrite the summations in ( |A.9|) , we obtain 



iVi 



D(e)de 
o exp(/3e) — 1 



C(3), T<T C , 



where 



C(3) 



00 x 2 dx 



1.202 • 



o e- 1 



(A.H) 



(A.12) 
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The transition temperature is defined by 



For T < T c , we have 



and therefore 



"=C(3)(^) 3 . (A.13) 



N = N + C(3) (^) , (A.14) 



^ = ^-(l)\ T<T C . (A.15) 



N \T, 



- c 



The exponent is changed from 3/2 (confined in a box) to 3 (confined in harmonic poten- 
tials). 

The total energy E in ( |A.6| ) is calculated in the thermodynamic limit [95|; iV — > oo, 
— > 0, iVu; 3 = fixed. The result is 

£ = 3WT^44, T>T C , 

J/sB ° 4 -C(4), T<T C , (A.16) 



(to) 



where with z = exp(/3/z) 



oo z l 



1=1 



9n(z) = Y,^ ( A - 17 ) 



and C(4) = vr 4 /90 = 1.082-. The specific heat C = dE/dT has a discontinuity, AC/Nk B = 
9((3)/((2) = 6.6. The transition is of second order. It is well-known that for a free boson 
in a box the derivative of the specific heat is discontinuous and the transition is of 3rd 
order (see for instance, Ref. [HI). 
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Appendix B 
Pseudopotential 



In this appendix, we give some details of deriving the pseudopotential. The argument 
follows Refs. [g [H§. 

We consider the two-body problem. Each particle has the mass, m, and an inter-particle 
potential, v (r), is the "hard-sphere" one, 

v(r) = < } ( (B.l) 

oo (r < a), v ' 

where a is the hard-sphere diameter with r the relative position vector between two par- 
ticles and r = \r\. The Schrodinger equation in the center-of-mass system is 

^(A + k 2 )tP(r)=v(r)^r), (B.2) 
where \i means the reduced mass, 

H = m/2. (B.3) 

Obviously, ip(r) is the wavefunction in the center-of-mass coordinate system, and (hk) 2 / (2/x) 
is the energy of the relative motion. Substituting ( |B.1| ) into ( |B.2| ), we have 

(A + k 2 )^)(r) = (r>o), 

ijj{r) = (r < a). (B.4) 

In terms of the spherical coordinate, 

r = (r sin 9 cos 0, r sin 6 sin <ft, r cos #) , (B.5) 
the solution of Eq. (|B.4|) for r > a can be written as 



oo +1 

*l>(r) = E E *im(0, 0) An (jiW - tan ^ n,(Ar)) , (B.6) 

1=0 m=-l 

with the boundary condition, 

?p(r)\ r=a = 0. (B.7) 
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Here Yi m (9, <fi) is a normalized spherical harmonic function, ji(x) and ni{x) the spherical 
Bessel and Neumann functions respectively, and A [m and r]i constants. We note that the 
constant rji is determined by the condition flB.7|) as 



tan 7ft = ji(ka)/rii(ka). 
The scattering length a/ for the partial Z-wave is defined by 



a t = —hmtanr]i(k)/k. 



(B.9) 



In what follows, we assume that the energy of the relative motion (hk) 2 /(2/x) is sufficiently 
small, and thus we consider a spherically symmetric (s-wave) solution, 



if)(r) = A(j Q (kr) - tanr/ n (kr)) 



where 



sinx/x, 
— cosx/x, 

A 00 /V47T. 



Jo{x) 
n (x) 
A 

From (|B~8|) , ( g3g ) and (gig ), we have 

tan 7/0 = — tan(/ca), 

leading to 

r/o = — ka. 



(B.10) 



(B.ll) 
(B.12) 
(B.13) 



(B.14) 
(B.15) 



Thus, a is identified with the s-wave scattering length. 

An idea of the pseudopotential is as follows: we find an equation with some "potential" 
such that (|B.10|) is the solution everywhere. For sufficiently small x, jo(x) and n${x) behave 
like 

j (x) w 1, n (x) w -l/x (x<l). (B.16) 
Thus, from (p.lOj ), for sufficiently small kr, we get 



which gives 



Remark that the relation 
jo(kr) satisfies 



r^r) =A(r+ ^) 



(B.17) 



(B.18) 



18|) is used only at r — 0. Because j'o( x ) is regular at x = 0, 
(A + k 2 )j (kr) = 0, (B.19) 
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for all r. On the other hand, n$(x) is singular at x — 0. Then, we calculate 

F (r) = (A + k 2 )n (kr), (B.20) 

with care. We integrate Fo(r) over a sphere V of radius e about the origin. From ( |B.20| ), 
we have 

/ d 3 rF (r)= / d 3 r An (kr) + k 2 [ d 3 rn (kr). (B.21) 
Jv Jv Jv 

By applying the divergence theorem to the first term in the right hand side of Eq. ( |B.21| ) , 

we get 



v 



d 3 r Ano(fcr) 



c)V 



dS ■ Vrio(kr) 







Aire -^n (kr) 

r=e 
47T 

47re sin(fce) + — cos(fce) 
k 



(B.22) 



The second term in Eq. (|B.21|) gives 



k 2 / d 3 rn (kr) = 4ir r 2 dr 



- cos(/cr)\ 
kr J 

— 47resin(fce) cos(fce) 

k 



Att 

T' 



(B.23) 



Substituting ( |B~22D and ( [B~23D into ( [B~2lD , we obtain 



47T 

d 3 rF (r) = — . 
v k 



(B.24) 



Noting that Fq(v) is identically equal to zero for r ^ 0, we conclude from (p.24j ) that 



47T 



F (r) = (A + k 2 )n (kr) = —5(r). 

k 



(B.25) 



Using (|B~T4|) , (pT8D , fpT9D and (|B~25| ) in Eq. ( plO|) , we have an equation that the solution 
([B.lOp satisfies everywhere, 



4tt 



d 



(A + k z )ip(r) = — tan(fca) 5(r) — (rip(r)) . 



k 



dr 



(B.26) 



For sufficiently small ka, we can replace tan(fca) by ka. Then, by dividing the both sides 
of Eq. (|B.26|) by /i 2 /(2/i), we finally arrive at 



h 



where 



^-Aif>(r)+v(r)ij>(r) = ^-k 2 ^(r), 
2fi 2/x 



Anh 2 a d 
v(r) = b(r) — (r ■) . 



dr 



(B.27) 



(B.28) 
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The operator v(r) ( |B.28| ) is known as the pseudopotential |5], [75| . We note that d/dr (r ■) 
appearing in (|B.28| ) is not a hermitian operator. But, if ip(r) is well behaved, namely 
differentiable at the origin, we can replace d/dr (r ■) by unity. So far, we have considered a 
to be positive. In general, however, the "diameter" of the hard-sphere a can be extended 
to be negative. This occurs when we may replace the low energy scattering from an 
attractive inter-particle potential of finite range by that from a hard-sphere one, known as 
the "shape-independent approximation." 



64 



Appendix C 

The Ground State Energy under the 
Thomas-Fermi Approximations 



We give some details of calculations of the ground state energy of the condensate under 
the one-, two- and three-dimensional Thomas-Fermi approximations. 

First, we consider the one- dimensional case. We start from the one-dimensional Ginzburg- 
Pitaevskii-Gross equation [69|-[T2"| with a harmonic potential term, 

h 2 d 2 ^ m 



Here and hereafter, g{> 0) means the strength of the inter-atomic interaction, and /i is 
the chemical potential. In the Thomas-Fermi approximation, the first term (the kinetic 
term) in the left hand side of Eq. (|C.1|) is neglected, which gives the number density of the 
condensate, |\1/| 2 , as 

|tff = I^_ ~mcu 2 x 2 y (C.2) 

By integrating the density (|C.2|) in the interval [— xo,xo\, where xo = (mu 2 ), we get 
the number of the particles N, 

N= f X ° dxm 2 = ^(^) 1/2 ^ 2 , (C.3) 



3g Vmw 



leading to 



Substituting ( |U.4j ) into the thermodynamic identity, 



*= W , (C5) 



we have the energy, E^d, given by 



£lD 3 m -^y% 5 , (a6) 
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Next, we consider the two-dimensional case. The two-dimensional equation is 

h 2 ( d 2 



2m 



(Q 2 Q 2 \ 771 



(C.7) 



According to the Thomas-Fermi approximation, we ignore the kinetic energy terms in the 
left hand side of (|C.7|) and get 

l^l 2 = - (a« - \muJ 2 x x 2 - \muj 2 y y 2 ^ . (C.8) 

By integrating ( |C.8| ) in the region /i — \mi>j 2 x x 2 — ^mu 2 y 2 > 0, we have the number of 
particles N, 

r%o ryo 

N= dx dy\^\ 2 , (C.9) 
J-x J —yo 

where x and y are defined as 



( 2w \ 1/2 ( 2 \ 1/2 / 1 \^ 2 
x = , y = - U - -mw 2 x x 2 . (C.10) 



Substituting ( Jfl% and ( |CT0l ) into ( |C^D , we have 



3g \mujy) \mui x j 
from which we get 

Using ( |C12| ) in ( |C5| ), we obtain the energy in the two-dimensional case, E 2 r>, 

E2D * fBy/V3,y/Y^y7^y%3 /2 (c . 13) 



3 V :>,;•,- / V 4 

Finally, we consider the three-dimensional case. The three-dimensional equation is 
Tt 2 ( d 2 d 2 d 2 \ Tn 

~ 2^ \e* + W + a? J * + 2 Ux + y + >* + m = (C 

As in the previous cases, we approximate the number density of the condensate as 

1/ 1 1 1 

9 



I^T = - - \™^ 2 x x 2 - ^muj 2 y y 2 - ^mu 2 z z 2 ^j . (C.15) 



The number of particles N is given by 



rxo ryo , 

N = I dx dy I dz\ty\ 2 , (C.16) 

xo J -yo 
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where xq, yo and z are defined as 



z = — » hu - tt™^ 2 - ^ muJ ly 2 ) ■ (C.17) 



Substituting ( |U.15|) and ( |C.17|) into ( |C.16| ), we obtain 



8vr 



1/2 / Q \ 1/2 / Q \ 1/2 



n= tt K2 hi h ^ (ai8) 

lop J \mujy J \muj£J 



which gives 

From ( |G.5| ) and (U.19), we have the energy in the three-dimensional case, denoted by E 3 n, 

To summarize, the ground state energy in the Thomas-Fermi approximation for the 
(^-dimensional Ginzburg-Pitaevskii-Gross equation with harmonic potential terms has the 
particle number dependence as 

£ rfD ~iV {(i+4)/(ci+2) . (C.21) 

This iV-dependence can be used to identify the effective dimensionality of the condensate 
under the anisotropic magnetic traps in the strongly repulsive case. 
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Table and Figure Captions 



Table Experimental values of parameters, a; s-wave scattering length, T c ; critical 

temperature, N t ; the total number of atoms in trap, lo\ trap frequency. Note that the total 
number of the atoms is the sum of normal and condensed ones. The values of T c , N t and 
u are taken from the reference in the right most column. 

Fig. |3.1.1| Approximations for the repulsive inter-atomic interaction: I. weak interaction 
case (|3.1.24] ) , II. strong interaction case ( |3.1.31| ), III. intermediate case-1 ( |3.1.40| ), IV. 
intermediate case-2 (|3.1.52| ). The abscissa and ordinate represent the anisotropy of the 
trap, 5, and the strength of the interaction, G±, respectively. 

Fig. |3.2.1| Stability of the two-component condensate: jVj = N 2 case for (a) a > 0, (b) 
a < 0. The boundaries of stable and unstable regions (the solid lines) are determined by 
Eqs. (|3.2.9[ ) and ( |3.2.10| ). The dashed lines are discussed in |3.2.4| . 

Fig. |3.2.2| Phase diagram for a n = a 22 = a > and a± 2 < —a. 

Fig. |3.2.3| Phase diagram for a u = a 22 = a < 0. (a) a 12 > 0, (b) a < a\ 2 < 0, (c) 

a.\i < a < 0. 

Fig. |3TX1| Function ^((3) defined in (|3.2.35|) . We set a n : a 22 : a X2 = 1 : -4 : 2 _ i. 

Fig. |3.2.5| limTVj^oo f{N\)/Ni as a function of 23 a^/o-u- We set an : a 22 = 1 : —4. 

Fig. |3.2.6| Phase diagram for an > and a 22 < 0. (a) a%2 > 2~^a; 11 , (b) < «i 2 < 
2~2a 11) ( c ) « 12 < 0. 

Fig. |4.2.1| N CjDew /N as a function of A; for the case that 5 = 0.718 [12 
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